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Abstract 

We propose a new systematic approach that allows one to derive 
the spin foam (state sum) model of a theory starting from the corre- 



> 
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sponding classical action functional. It can be applied to any theory 
whose action can be written as that of the BF theory plus a func- 
tional of the B field. Examples of such theories include BF theories 
with or without cosmological term, Yang-Mills theories and gravity in 
various spacetime dimensions. Our main idea is two-fold. First, we 
propose to take into account in the path integral certain distributional 
configurations of the B field in which it is concentrated along lower 
dimensional hypersurfaces in spacetime. Second, using the notion of 
generating functional we develop perturbation expansion techniques, 
with the role of the free theory played by the BF theory. We test our 
approach on various theories for which the corresponding spin foam 
(state sum) models are known. We find that it exactly reproduces the 
known models for BF and 2D Yang-Mills theories. For the BF theory 
with cosmological term in 3 and 4 dimensions we calculate the terms 
of the transition amplitude that are of the first order in the cosmologi- 
cal constant, and find an agreement with the corresponding first order 
terms of the known state sum models. We discuss implications of our 
results for existing quantum gravity models. 



e-print archive: [http://xxx.lanl.gov/abs/hcp-th/9807092 
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1 Introduction 



Recently Baez jlj] , motivated by results [|j of the canonical (loop) approach to 
quantum gravity, introduced a notion of "spin foam models" . The spin foam 
model of a theory, as we explain in details below, gives a way to calculate 
transition amplitudes for this theory. A transition amplitude arises as a sum 
over branched colored surfaces, or, using the name proposed by Baez, spin 
foams. The results of Reisenberger and Rovelli Reisenberger ||, Barrett 
and Crane |||] indicate that quantum gravity might be constructed as a spin 
foam model. 

However, spin foam models arise not only in the context of canonical 
(loop) quantum gravity. Indeed, many of the known TQFT's — for example, 
the Migdal-Witten model in two dimensions ^] , the Ponzano-Regge R and 
Turaev-Viro ||] models in three spacetime dimensions, the Ooguri model § 
and its generalizations by Crane, Kauffmann and Yetter [10] — are state sum 



models and, by considering instead of sums over colored triangulations sums 
over the corresponding colored dual complexes, can be formulated as spin 
foam models. The state sum models mentioned are believed to correspond 
to certain well-known topological field theories. However, a precise relation 
between the classical theory and the corresponding state sum model is not 
known for some of these models. Moreover, even for the models whose 
relation to the corresponding classical theory is established, this relation is 
rather indirect. Thus, summarizing, one can say that there exists a large 
class of spin foam models that are claimed to correspond to such classical 
theories as Yang-Mills, BF theories and gravity in various dimensions, but 
no precise relation to this effect is known in many cases. 

The present paper is devoted to the study of this relation. More precisely, 
we study the relation between a spin foam formulation of a theory and the 
corresponding classical action principle. We propose and study a particular 
procedure that allows one to get a spin foam formulation for a large class 
of theories starting from the corresponding action principle. Our approach 
is not limited to any particular spacetime dimension. However, having in 
mind applications to the models mentioned above, we restrict ourselves to 
two, three and four spacetime dimensions .[[] 

When considering gravity, for simplicity, we restrict ourselves to metrics 
of Euclidean signature. 

Our approach is inspired first of all by the results of the canonical (loop) 



approach 112, 13, 14, lq] to quantum gravity. Of great importance for us were 



also the ideas of Regge [16], the work of Baez [17] on degenerate solutions 



of general relativity and some results of Barrett [18]. Motivated by all these 



1 For an application of our procedure to higher- dimensional gravity see |jnj. 
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results we propose to "approximate" some of the continuous fields appearing 
in the Lagrangian of the theory by distributional fields concentrated along 
lower dimensional hypersurfaces in spacetime. To make contact with the 
known lattice models, we put our distributional fields on a triangulated 
manifold, so that they are concentrated along polygons of the dual complex. 
However, the idea is also valid in a context more general than a lattice model. 
The action of the theory calculated on such distributional fields becomes a 
sum of certain integrals over the hypersurfaces of lower dimension. Then, 
to calculate a transition amplitude for the theory, one constructs the path 
integral of the exponentiated action over such distributional fields. As we 
shall see, this path integral becomes quite tractable. 

To make this idea little more precise, let us specify the class of theories 
that can be treated within our approach. We consider theories whose action 
is BF-like, that is, whose action is of the form 



where E is the Lie algebra valued two-form field (for instance, B field of BF 
theory), F is the curvature of the connection form and $ is certain (poly- 
nomial) function of the E field, which can also depend on some Lagrange 
multipliers, as in the case of gravity in four dimensions, or on an additional 
background structure, as, for example, a fixed measure on M. in the case of 
Yang- Mills theory in two dimensions. Below we will give many examples of 
theories belonging to this class. Thus, the action is that of the BF theory 
with an additional term. We call the B field E because of its relation with 
the non-abelian "electric" field of the canonical formulation. 

In order to compute a transition amplitude for a theory belonging to 
this class, we take the path integral over classical configurations where the 
E field, which is a zero-, one- or two-form in different dimensions, is a distri- 
butional field that vanishes everywhere except on two dimensional polygons 
in spacetime. Note that in our approach only these distributional configu- 
rations of the E field are taken into account in the path integral. No usual 
smooth configurations of the E field are summed over. This can be viewed 
as a discretization procedure where we "approximate" the smooth E field by 
"squeezing" it into a distributional field on a collection of two-dimensional 
polygons. One can expect that this approximation becomes good as the two- 
dimensional polygons become sufficiently packed in spacetime. However, this 
procedure can also be justified at a deeper level, as we shall argue. 

First, a justification for considering distributional fields in the path in- 
tegral comes from the canonical quantum theory. The phase space of all 
theories we work with, as it can be seen from the above action functional, 
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is that of Yang-Mills theory. The space of quantum states of the canonical 
theory is then given by L 2 (A/Q), where A/Q is the space of connections on 
the spatial manifold modulo gauge transformations, and L? is defined by an 
appropriate choice of measure on this space. In the case the spatial manifold 
is one-dimensional, there exists quite a natural choice of this measure: the 
functionals on A/Q are just class functions, and it is natural to use the Haar 
measure on the group to define the inner product of such functionals. In 
higher dimensional cases, an analogous construction of L 2 (A/G) was devel- 



oped [14] within the approach of canonical (loop) quantum gravity. In this 
construction, integrable functions arise by considering the so-called cylindri- 
cal functionals, that is, functionals that depend on connection only through 
holonomies along some paths in space. In other words, one constructs the 
space of states of the theory by considering spaces of states of lattice gauge 
theories on graphs in space, and then taking the projective limit over graph 
refinements (see [14] for details). A characteristic feature of the quantum 



theory constructed with such a space of states is the distributional nature 
of the electric field operator. Let us illustrate this on the case of (2+1) di- 
mensional theory. A typical state from L 2 (A/G) depends on the connection 
only through holonomies along some paths in space. Although this is not 
important for our heuristic consideration here, let us note that gauge invari- 
ance implies that the paths form a closed graph in space. Let us denote the 
edges of such a graph by e«. Then a typical state supported on this graph is 
of the form 

* (h ei (A), . . . , h en (A)) , 

where h ei (A) is the holonomy of A along the edge ej, and \P is a function 
with sufficiently good behavior. The electric field E, which classically is 
a canonically conjugate quantity to A, quantum mechanically becomes the 
operator 5/5A(x) of variational derivative with respect to A. This operator, 
when acting on a typical state, gives zero for all points x except the points 
lying on the edges e^, where the result is distributional. Thus, in a typical 
quantum state, the electric field E is distributional with support on edges 
ej of the corresponding graph. One normally views a state of the canonical 
theory as a state "at a given time". Thus, electric field "at each given time" 
is distributional and concentrated along edges of the graph. This suggests 
that one should think of the "history" of the electric field configuration as of 
a collection of two-dimensional surfaces in spacetime, along which the elec- 
tric field is distributional. When sliced by a spatial hypersurface "at given 
time" , such a collection of surfaces gives a graph in space, and electric field is 
distributional along the edges of this graph. These heuristic considerations 
suggest that one should consider configurations of E field that are distribu- 
tional and concentrated along two dimensional surfaces in spacetime. As we 
will see in more details below, the surfaces of support of E field are always 
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Figure 1: A typical quantum state of (2+1) dimensional theory is labelled 
by a graph in space. The corresponding functional of connection depends on 
A only through holonomies along edges ej. The electric field in such state is 
distributional and concentrated along the edges e^. 



two-dimensional, independently of the dimension of the spacetime manifold. 

The second justification for considering the path integral over such dis- 
tributional E fields comes from the fact of existence of similar distributional 
solutions of the equations of motion in the classical theory. As it was shown 
by Baez [|l7j for the case of gravity, there exist solutions of the classical 
equations of motion in which the E field is distributional and concentrated 
along two-dimensional surfaces in spacetime. In these solutions the E field 
is zero except in a tubular neighborhood of a two-dimensional surface in 
spacetime. The connection field away from the surface then must be either 
arbitrary or flat, depending on the detailed form of the Lagrangian of the 
theory. Solutions of this type can be shown to exist for a large class of theo- 
ries. Here we take the fact of the existence of such solutions as a motivation 
for our construction. Namely, we know that the simplest way to approach 
the computation of a path integral is to consider the semi-classical approx- 
imation. Then the Feynman path integral is evaluated as an integral over 
the classical moduli space (space of classical solution) with a measure given 
by the one loop expansion of the action around the classical solution. In this 
case, there may be different phases in the theory, and the type of solutions 
that dominate the path integral may be different in different phases. Thus, 
there may exist a phase in a theory of the type we study where the path 
integral is dominated by the distributional solutions found by Baez. At this 
point it is instructive to make an analogy with the usual Yang-Mills theories. 
In (abelian) Yang-Mills theories one also finds "distributional" solutions of 
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the classical equations of motion: monopole solutions. In Dirac's way [|19| of 
treating the monopole, it is described by a configuration of fields in which 
magnetic field is distributional along certain lines. Distributional monopole 
solutions are not in the spectrum of perturbative quantum theory. However, 
one must include them when calculating the partition function of the theory 
via path integral method to reproduce correctly properties of certain phases 
of the theory, such as the confining phase. In fact, the confining phase of 
the theory is the one in which monopole solutions dominate the path inte- 
gral |2(J. Distributional solutions of the type discovered by Baez [|l7| are 
very similar in nature to monopole solutions. In fact, they are close to dual 
monopoles: electric field is distributional and concentrated along flux lines 
in space, instead of the magnetic field. Thus, the analogy with monopole so- 
lutions suggests that one should include these distributional solutions when 
taking the path integral. In this paper we include only these distributional 
solutions, neglecting the usual smooth solutions playing the role in pertur- 
bative quantum theory. Thus, we attempt to study the phase of various 
theories in which the distributional solutions dominate the path integral. 
This is the motivation for considering the path integral over distributional 
configurations of E field. 

The last but not least justification of our approach is the test of the 
procedure on the exact results coming from TQFT's. As we shall see, for a 
large class of theories our approach gives results that are in a good agreement 
with the known state sum models. 

In this paper we only calculate path integrals over the distributional 
fields "living" on a fixed collection of surfaces in spacetime. In the case of 
topological theories, the result for a particular fixed collection of surfaces 
does not depend on the collection. One cannot expect this property to hold 
for such theories as gravity, which is not a topological theory. Thus, one 
might want to perform a sum over all possible collections of surfaces in 
spacetime, or take a limit as the triangulation becomes more refined. We do 
not address this, very important, problem in the present paper. 

Usual interactive quantum field theories have been successfully under- 
stood using perturbative expansion in terms of Feynman graphs. In this 
case the interacting QFT is considered as a perturbation of a free field the- 
ory. From a technical point of view, what makes this approach successful is 
the exact solvability of the free field theory, i.e., our ability to compute all 
possible correlation functions. This ability can be recasted into the knowl- 
edge of the generating functional, and Feynman graphs can be viewed as the 
evaluation of certain differential operators acting on the generating func- 
tional. 
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We will employ a somewhat similar strategy. In order to calculate tran- 
sition amplitudes for our class of theories we will use the machinery of gen- 
erating functionals. In our approach the role of the "free" field theory is 
played by the BF theory and the terms in the action polynomial in E are 
analogous to the "interaction" terms. The idea to use the BF theory as a 
"free" field theory is not new: it has been applied with some success in the 
context of Yang-mills theories (see f2lj). However, the details of how the BF 
theory is used in our approach differ from those of 21]. In order to calculate 



J=o 



a transition amplitude, i.e., the path integral of the exponentiated action 

J VAVE e l / IMEaf)+^(E)}^ 

we formally rewrite this path integral as 
where 

Z[J] := J VAVE e ifmEAF)+Tr(EAJ))_ 

We shall refer to Z[J] as the generating functional. Here J is a Lie algebra 
valued two form field. Using terminology from field theory we will call the 
J field current. One of our main results is the exact computation, in the 
context of spin foam models, of the generating functional Z. It is obtained 
by integrating over fields A, E that "live" on a fixed triangulation of the 
spacetime manifold. Then the transition amplitude for any theory of the 
type we consider is given as a formal power series in variational derivatives 
5/5 J with respect to the current field. The series we get are quite reminiscent 
of the usual Feynman diagram expansions. Thus, the powerful technique of 
generating functionals allows us to study different theories — such as BF 
theory and gravity — from the same point of view. We note that, to the 
best of our knowledge, the technique of generating functionals has never 
before been used in the context of spin foam models. 

In cases when the spacetime manifold is two, three and four-dimensional, 
for each fixed triangulation we find an explicit expression for the generating 
functional. We then use the obtained expressions to calculate the transi- 
tion amplitudes for various theories. We find that our approach exactly 
reproduces the state sum models of 2D Yang-Mills and BF theory in any di- 
mension. For BF theories with cosmological term we compare the amplitude 
obtained within our approach with the one given by the corresponding state 
sum model by decomposing them into the power series in the cosmological 
constant and comparing the first order terms. Again, we find an agreement. 
We make some comments as to consequences of our results for the quantum 
gravity models in four spacetime dimensions. 
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The paper consists of two main parts. In the first part, containing this 
section and Sec.||, we introduce the path integral techniques in the context 
of spin foam models. Here we also introduce the generating functional, 
and compute it in various spacetime dimensions. The second part contains 
applications of the general techniques developed earlier to concrete theories. 
This also serves as a test of our approach, for one can compare the spin 
foam models obtained by our techniques with the known state sum models. 
In Sec.|| we define the class of theories that can be treated by techniques 
described in this paper, and give examples of physically interesting theories 
belonging to this class. Sec.^| describes briefly some known state sum models 
and summarizes the known arguments that relate the state sum models to the 
classical theories. Here we also obtain some properties of various state sum 
models that will be needed later. Finally, using the generating functional 
technique, we derive in Sec. | the spin foam models corresponding to various 
theories of Sec.|| and compare them with the known models. We conclude 
with a discussion. 



2 Generating functional 

As we briefly explained in the Introduction, to calculate transition ampli- 
tudes for a theory of the type we consider here it is very convenient to first 
calculate the generating functional of BF theory. Then transition ampli- 
tudes can be obtained as formal power series in derivatives with respect to 
the current. In this section we calculate the generating functional of the BF 
theory, i.e., the path integral 

Z\J\ = J VAVE e ifmEAF)+Tr(EAJ)] ^ 

in two, three and four spacetime dimensions. Here A is a connection on 
some G bundle over the spacetime manifold Ai, where G is the gauge group 
of the theory, F is the curvature 2-form of the connection, and E, J are 
Lie algebra valued (d-2)-form fields, where d is the dimension of Ai. For 
simplicity, we consider here only the case when the spacetime manifold Ai 
is closed. We also restrict ourself to the case of the gauge group being 
G = SU(2). Although for practical applications of spin foam models one is 
usually interested in the case when Ai has a boundary consisting of "initial" 
and "final" spatial hypersurfaces, we do not consider this here because our 
main goal is to compare the spin foam models obtained by our techniques 
with the known state sum models, which can be done already in the case of 
closed Ai. At any rate, the inclusion of boundaries is straightforward, and 
all our formulas can be easily generalized to this case. To give meaning to the 
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path integration over spacetime fields A, E, we replace the E field by certain 
distributional field concentrated over two-dimensional surfaces in spacetime. 
To make contact with the known state sum models we put our fields on the 
2-dimensional cellular complex dual to some triangulation of the spacetime 
manifold. In each particular calculation this triangulation is fixed and all 
results explicitly depend on it. We first investigate the generating functional 
in a general case, without specifying the dimension of spacetime manifold, 
and then find an explicit expression for it in each particular dimension. 

2.1 General framework 

Let us fix a triangulation A of a d-dimensional compact oriented spacetime 
manifold M. This triangulation defines another decomposition of M into 
cells called dual complex. There is one-to-one correspondence between k- 
simplices of the triangulation and (d-k)-cells in the dual complex. We orient 
each cell of the dual complex in an arbitrary fashion, which also defines an 
orientation for all simplices of the triangulation. The (d-2)-form E can now 
be integrated over the (d-2)-simplices of the triangulation the result being a 
collection of the Lie algebra elements X — one vector X € su(2) for each 2- 
cell dual to a (d-2)-simplex of the triangulation. We would like to discretize 
our model by replacing the continuous E field by the collection of the Lie 
algebra elements X. It turns out, however, that this is not yet the most 
convenient set of variables for the theory. For reasons which we give below, 
we use another, more convenient set of variables introduced by Reisenberger 

In what follows we will often call 2-cells of the dual triangulation dual 
faces. Instead of a single X for each dual face, we introduce a set of vari- 
ables, which we call X w . To do this, we divide each dual face into what 
Reisenberger calls "wedges" (see Fig. 2). To construct wedges of the dual 
face one first has to find the "center" of the dual face. This is the point 
on the dual face where it intersects with the corresponding (d-2)-simplex of 
the triangulation. One then has to draw lines connecting this center with 
the centers of the neighboring dual faces. The part of the dual face that lies 
between two such lines is exactly the wedge. Thus, each dual face splits into 
wedges, and we assign a Lie algebra element X w to each wedge w. All vari- 
ables X w that correspond to a single dual face are required to have the same 
length. Recall that vectors of the same length only differ by a gauge trans- 
formation. Wedges of a given dual face are in one-to-one correspondence 
with the d-simplices of the triangulation neighboring this dual face. Thus, 
the physical meaning of each variable X w can be said to be the integral of 
E over the (d-2) simplex of the triangulation "from the point of view" of a 
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Figure 2: A face a of the dual triangulation. The portion of a indicated by 
bold lines is what is called wedge here. The point labelled by C is the center 
of one of the d-simplices neighbored by a. 

particular d-simplex containing this (d-2)-simplex. In other words, with the 
introduction of the wedge variables X w we allow integrals of E over one and 
the same (d-2)-simplex of triangulation to be different when this simplex is 
considered as belonging to different d-simplices, as long as lengths of all such 
X w coincide. Note that the number of variables X w that arises this way for 
a given triangulation is equal to the number of d-simplices times the number 
of (d-2)-simplices in each d-simplex. 

Having discussed the geometrical meaning of the wedge variables X w 
we are ready to introduce the distributional E fields. Heuristically, our 
procedure of replacement of a smooth E field by a distribution concentrated 
along the wedges amounts to "squeezing" of the smooth E field that is 
"spread" over a (d-2) simplex of the triangulation to a single point on this 
(d-2) simplex, the point where the simplex intersects with the 2-cell of the 
dual complex. Thus, we define a distributional field E w concentrated along 
a wedge w to be a 2-form satisfying the following relation: 



Here J is any Ad(P)-valued 2-form, and w stands for a wedge. The integral 
on the right-hand-side is performed over the wedge w. This, in particular, 
implies that the Lie algebra element X w is equal to the integral of E w over 
the (d-2)-simplex of the triangulation that is dual to the 2-cell a. We then 
define distributional field E to be 






(3) 



w 



Although with the introduction of the wedge variables we have increased 
the number of independent variables in the theory, one can still argue that 
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the physical content of the theory should be the same because the norm of 
all the X w belonging to one and the same dual face are required to be equal. 
For the case of BF theory one can directly prove that the introduction of 
the independent wedge variables does not change the theory in the case the 
topology of dual faces a is that of a disc, which is what we will assume here. 
For more complicated theories with non-zero "interaction" term (the term 
in the action that makes the theory to be different from the BF theory), 
the wedge variables are so indispensable that it is not clear whether one can 
even define the theory without these variables. It is interesting, however, 
that, even in the case of BF theory one has to use the wedge variables if 
one wants to reproduce a triangulation-independent state sum model for the 
case the topology of a is different from that of a disc. This serves as a 
strong motivation for considering the wedge variables. Another motivation 
for using them is, as we explain below, that the wedge variables appear 
naturally when one tries to implement the gauge symmetry in the theory. 
Also, as we shall see, using these variables we will be able to reproduce the 
well-known state sum models in various dimensions. Thus, we are going to 
use the set of wedge variables motivated, in particular, by the results we 
obtain. 

The discretization we consider is consistent with the results of the canon- 
ical approach. The restriction of E constructed as described above to a 
generic "spatial" ((d-l)-dimensional) slice E of M. is non-zero only along 
edges of the graph r in E, where T is the intersection of the dual complex 
with E. Thus, this "discretization" realizes rather explicitly Faraday's idea 
of lines of electric force. Indeed, the restriction of E on a spatial hyper- 
surface can be described as a collection of flux lines of non-abelian electric 
field. It is not hard to see what states would correspond to such E in the 
quantum theory. Recall that the quantum states are functionals ^ [A] , where 
A is the pullback of the connection field A on the spatial hypersurface E. 
The densitized vector field E (that is dual to E) becomes the operator of 
variational derivative with respect to A. In order for E in a state ^[A\ to 
be zero everywhere except along the edges of T the state ^ should depend 
on the values of the connection A only along the edges. Note also that our 
construction of the distributional E field gives us E whose length is constant 
along each edge. This means that the state $ can depend on the values of 
A only as a function of the holonomies of the connection along the edges. 
These are exactly the states one finds in the canonical approach to quantum 
theory, as we explained in the introduction. 

To calculate the generating functional Z as a function of a fixed trian- 
gulation A, we have to take the integral of the exponentiated action over 
the "discretized" dynamical fields A, E, that is, the fields "living" on the 
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Figure 3: A wedge w of the dual face and the group elements: holonomies 
along the edges of w. 

triangulation A. The "discretized" action, that is, the action evaluated on 
the distributional field E, becomes a function of the Lie algebra elements X, 
and a functional of connection A and current J. Using (||),(|3|), we find that 
this discretized action is given by 

]T / [Tt(FX w ) + Tt(JX w )}, (4) 
w Jw 

where the sum is taken over the wedges of faces of the dual complex ('w' 
stands here for a wedge), and the integral is performed over each wedge, the 
integrand being the curvature of the connection A contracted with the Lie 
algebra element X w "living" on that wedge plus the current J contracted 
with X w . Each integral is performed using the orientation of the dual face 
to which the wedge belongs. 

We could now substitute this discretized action into the path integral and 
integrate over X and A. However, we first have to discuss what measure has 
to be used to integrate over X, A. To give meaning to the integration over 
A, let us replace a continuous field A by a collection of group elements. To 
do this we use the following approximation 

/ Tr(F X w ) w Tr(Z w X w ), (5) 

where Z w is the Lie algebra element corresponding to the holonomy of A 
around the wedge (see Fig. 3) . The base point of the holonomy is not fixed 
at this stage (see below). In other words, 

exp Z w = g 1 h 1 h 2 g2, (6) 

where gi,hi,h,2, Q2 are the holonomies of A along the four edges that form 
the boundary of the wedge w. We assume a local trivialization of the bundle 
over w so that the holonomies are group elements. The order in which the 
product of group elements is taken is determined by the orientation of the 
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dual face. Of course, exp Z w is defined only up to its conjugacy class (any of 
the four edges can be taken to be the "first"). Thus, there is an ambiguity in 
the choice of the base point for the holonomy, which we have to fix in some 
way. To fix it we define a notion of "discretized" gauge transformation and 
fix the ambiguity requiring the discrete action to be gauge invariant. 

But first, let us replace the continuous current field J by a collection of 
Lie algebra elements J w : 

J w := I J- (7) 

J w 

Here, to perform the integration, a trivialization of the bundle over the wedge 
w is chosen. The discretized action now becomes 

J2^(Z w X w ) + Tr(J w X w ). (8) 

w 

We fix a definition of Z w in such a way that this action is gauge invariant. 

We define the discrete gauge transformation so that it "acts" in the center 
of each d-simplex. More precisely, a gauge transformation is parameterized 
by a collection of group elements: a group element g for each d-simplex. 
First of all, the gauge transformation is defined for the holonomy U of the 
connection along any loop that starts and finishes at the center C of the 
d-simplex. The transformation is as follows: 

U^gUg- 1 . (9) 

The wedge variables X w and the discrete current variables J w transform as 

X w -> gX w g~ l , J w -> gJ w g v . (10) 

With this definition Tr(X w Z w ) is gauge invariant only when exp Z is defined 
as the holonomy around the wedge w whose starting and final point is the 
center C of the d-simplex, as in (||). This fixes the ambiguity in Z w . With 
Z w defined this way the discretized action is gauge invariant. 

The approximation (||) is good for Z w close to zero element in the Lie 
algebra. Thus, this approximation is certainly justified for BF theory, where 
we expect only connections close to flat to matter in the quantum theory. 
It is harder to justify this approximation for theories, for which the classical 
equations of motion do not imply the connection to be flat, as is the case, for 
example, for BF theory with cosmological term or for gravity. However, even 
for such theories, one would expect the approximation (||) become better as 
the triangulation of the manifold becomes finer and the dual faces become 
smaller. In quantum gravity we expect finer triangulations to matter most, 
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which can serve as a justification for the above approximation in the case 
of gravity. The approximation (||) will be additionally justified when we 
discuss the problem of the integration over X w variables. To define the later 
we will refer to a certain standard field theory calculation for BF theory in 
two dimensions. As we shall see, the approximation (||) is quite natural from 
the point of view of 2D BF theory. 

The approximation (||) finishes the discretization procedure for the clas- 
sical action. We can now calculate the path integral for the generating 
functional Z by integrating the exponentiated discrete action (|8|) over the 
Lie algebra elements X w and group elements g, h. This path integral is given 
by 

Z(J, A) = TT / dgT] ( dhfjldXve* £^P«+^-M 

JSV(2) J -- L JSU(2) J V 

(11) 

Here the integrals are taken over all group elements g, h, entering the expres- 
sion through Z w , see (0). These integrals form the discrete analog of DA, 
dg is the normalized Haar measure / dg = 1 on SU(2). The integrals over 
X w present here - one for each edge - form the analog of the integral over 
DE. The measure dX here is some measure on the Lie algebra. For now, 
we will leave this measure unspecified. Recall now that the integrals over 
X w are not independent: all X w that belong to one and the same dual face 
should have the same norm. However, as we shall see, there is no need to 
impose these constraints: one can integrate over all X w independently, but 
later, when one integrates over the group elements, only the contributions 
that come from X w of an equal norm will survive. Thus, we do not impose 
these constraints in ([Tl]) . 



Let us now investigate the structure of the path integral (|l l[) . To calcu- 
late Z we have to find the function of exp Z w , exp J w that is given by 

dX w e iTv ( XwZw+XwJw \ (12) 

In fact, it is not hard to see that this function is proportional to the (^-function 
of exp Z w peaked at exp J w . The proportionality coefficient can be a gauge 
invariant function of J w . As we explain below, this function must be set to 
be equal to P{J W ), where P is the function that relates Lebesgue measure 
on the Lie algebra and Haar measure on the group (see the Appendix B). 
Thus, as the result of (|l^) we get: 

P(J w )8(expZ w expJ w ), (13) 



where 5 is the standard ^-function on the group. Thus, the function of 
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Z„, J w given by ( |12| ) can be written as 

P(J w )^2dimjXj(eyLpZ w expJ w ), (14) 

3 

where we used the well-known decomposition of 5-function on the group 
into the sum over characters Xi(exp Z w exp J w ) of irreducible representation. 
Here the sum is taken over all irreducible representations of SU(2) (labelled 
by spins j), and dimj = (2j + 1). 

Thus, we find that 
z ( J ' A ) = n ^11 / dh Y[ p ( J w)J2 dim 3\nXj w ( e wZ w exp J w ). 

1L JSV(2) 1L JSV(2) yj~ ^ 

(15) 

The integration over the group elements can now be easily performed using 
the well-known formulas for the integrals of the products of matrix elements. 
However, let us first give a systematic derivation of the result (|l^). 



To find the result of the integral (|12|) we relate it to a more complicated 
integral, for which a precise result is known from 2D BF theory calculation. 
Recall that so far the measure dX w in (|l^) is unspecified. The following 
relation to 2D BF theory will also serve the purpose of specifying this mea- 
sure. Let us restrict our attention to a particular wedge w. We can restrict 
the bundle V to w and get an SU(2) bundle V w over w. Let A be a connec- 
tion on V w , and E be an Ad("P) valued 0-form. Consider the following path 
integral: 

J VAVEexpfi j [Tt(F E) + Tr(J£)]^ , (16) 

where the integration over A is performed subject to the condition that the 
connection on the boundary of w is fixed, and J is given by 

J = 5(p)J w , (17) 

where p is an arbitrary fixed point on w, and J w is the same as in (|l2]) , (p^) . 



The path integral (16) is just a partition function of BF theory on the disk 
with the distributional source given by (|i"7|). This partition function can be 
derived using results of |23| . The result is given by ( |i~3| ) , where Z w is the Lie 
algebra element that corresponds to the holonomy of A along the boundary 
of w. We will not present this calculation here. Instead, we refer to the 
calculation performed in 23] for the partition function of 2D BF theory on 



a punctured sphere. The result ( |13|) can then be checked by taking the 
partition function on the disk (equal to 5(g), where g is the holonomy along 
the boundary of the disk), and integrating it with (^) over dg. This must 
give the partition function on a punctured sphere, and indeed reproduces the 
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result given in . The only cautionary remark we have to make is that the 
calculation performed in [^] finds a gauge invariant partition function, that 
is, the one in which one takes J = 5(p)hJ w h^ 1 and integrates over dh. The 
techniques developed in pj| can be used only to calculate gauge invariant 
quantities, and are not directly applicable to the integral (plf). Thus, strictly 
speaking, using the results of |23| one can only argue that fll6D is equal to 

P{J w )5{e Zw he Jw h- 1 ), (18) 

where h is some group element. To get rid of h in this expression, and, thus, 
to get (|13|), we will recall how "discretized" gauge transformations act on 
the Lie algebra elements Z W ,J W , see (p|) , (|To|) . The result of (|l2|) must be 
invariant under this gauge transformations. It is not hard to see that this 
fixes h above to be unity, thus, giving (fH|). 

Having discussed how one can calculate the path integral (|l6|), let us 
now show that this path integral is, in fact, equivalent to (|T2|). Indeed, 
the integration over E in ( |l~6|) can be performed in two steps. First, one 
integrated over E(x),x ^ p, then one integrates over E{p): 

f dX w [ VE f VAexp (i f [Ti(F E) + Ti(J E)]) . (19) 

J JE(p)=X w J \ Jw J 

Since the current is distributional and concentrated at point p, the last term 
in the exponential does not matter when one integrates over VE,VA. On 
the other hand, using the same approximation as the one used in (||), it is 
not hard to show that 

/ VE f VAexp (i [ Tr(FE)) , (20) 

JE(p)=X w J \ Jw J 

where the integral over A is taken with A on the boundary of w fixed, is 
approximated by 

where exp Z w is the holonomy of A along the boundary of w. Putting this 
back to (|l^) one gets exactly (|l^) . This finishes our discussion of the deriva- 
tion of (111). 



Having discussed the derivation of the expression (15) for the generating 
functional, we can now perform the integrals over the group elements. In- 
tegration over the group elements h that correspond to edges dividing dual 
faces into wedges is the same in any dimension, and we perform it here. The 
rest of the group elements corresponds to edges that form the boundary of 
the dual faces. The integration over these group elements g is different in 
different dimensions, and we perform it in the following subsections. Each 
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Figure 4: The product of the group elements is the holonomy around the 
dual face a with the insertions of the group elements exp J at each center of 
the corresponding d-simplex. 

group element h is "shared" by two wedges; thus, we have to take the in- 
tegral of the product of two matrix elements. Such an integral is given by 
( |142| ) of the Appendix. Integrating over all these edges, and making a simple 
change of variables to eliminate trivial integrations, we find 



Z(J,A) = n/ dg e x (22) 

JSV(2) 



/SU(2) 

]J^(dim i J KCT P(Ji) • • -P(Jn)Xj CT (exp J x g tl exp J 2 ■ • - exp J n ge n )- 

Here the remaining integrals are over the group elements g e that correspond 
to the edges e of the dual complex (edges that connect centers of 4-simplices). 
The second product is taken over the dual faces a; j a is the spin labelling 
the dual face <r, K a is the Euler characteristics of a. It is equal to unity 
in case the dual face has the topology of a disc. In what follows we will 
always assume that this is the case. The order of the group elements in the 
argument of Xja ls clear from the Figure 4. 

Note that in the integration we just performed survived only the terms 
for which the spins of all wedges corresponding to one and the same dual face 
are equal. In (p2j) this spin, which is the same for all wedges belonging to 



the same dual face, is denoted by j a . As one can see from the formula ( 137|) 
for the inverse Fourier transform, the spin labelling a wedge in (|l^) has the 
meaning of the length of the corresponding Lie algebra element X. Thus, as 
we said above, there is no need to impose the constraint that length of all 
wedge variables X w belonging to one and the same dual face are equal: this 
constraint gets imposed automatically when the integration over the group 
elements is performed. 



The expression fl22|) for the generating functional can be further sim- 
plified by integrating over the group elements g e . However, the result of 
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this integration is different in different dimensions, and we will perform it 
separately for each particular case. 

Before we calculate the generating functional for each particular space- 
time dimension, let us pause to discuss some general properties of the ex- 
pression for Z(J, A). The generating functional has several symmetries 
which we describe as follows. First, the generating functional is invari- 
ant under gauge transformations. Namely, let us consider a transformation 
Jw 9h(w)Jw9i l ^ w y where h(w) denotes the d-simplex of A to which w be- 
longs. Here gh( w ) ls the same for all w belonging to the simplex h. The 
generating functional satisfies: 

Z(gJg-\A) = Z(J,A). (23) 



This "discrete" gauge transformation is parameterized by one group element 



for each d-simplex of A. However the expression (22) for the generating 
functional has a bigger invariance. Namely, let us denote by e+(w),e~(w) 
the two edges of the wedge w which meet at the center of the d-simplex h(w). 
Let us associate with each edge e±(w) a group element g e± ( w ) so that it is one 
and the same for different wedges w when e±(w) is one and the same. The 
transformation of all currents according to J w — > g e+ ( w )J w g~ 1 ^ leaves the 
generating functional invariant. This transformation, which is parameterized 
by d + 1 group elements per simplex, contains the gauge transformation as 
its particular case. The later corresponds to g e+ i w ) = g e _( w ) = 9h(w) f° r an 
w that belong to h. The appearance and significance of the described extra 
symmetry is not yet clear to us. 

The generating functional is covariant under the diffeomorphism group. 
Let / denote a diffeomorphism of the underlying manifold M , and denote by 
/(A) the image of the embedded triangulation A under the diffeomorphism. 
If J w = I w J, let us denote f*J = J f{w) J = J w f*J. Then 

Z(J,f(A)) = Z(f*J,A). (24) 



The generating functional is generally not invariant under a refinement of 
the triangulation A unless the theory is a topological field theory (unless 
J = 0). 

To define the generating functional we had to choose an arbitrary orien- 
tation of the wedges. However, the generating functional is independent of 
the orientation chosen. If one, for example, reverses the orientation of the 
wedge w this results: (i) in the change J w = J w J into — J w ; (ii) it changes 
the holonomy of the connection along the wedge into its inverse. The two 
effects cancel each other leaving the generating functional invariant. 
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Figure 5: Part of a triangulation of 2-dimensional spacetime manifold, with 
one face of the dual complex shown. Bold lines indicate a wedge of the dual 
face. 



2.2 2 dimensions 



The case of two spacetime dimensions is somewhat special in the sense that 
dual faces (2-cells) cover the manifold (see Fig. 5). The general construction 
of the previous subsection prescribes to replace the continuous field E, which 
in the case of two dimensions is a zero form, by a distributional field concen- 
trated along dual faces. However, since the spacetime is two-dimensional, 
the dual faces cover the manifold, and we just have to replace the E field by 
a field constant on each wedge, and equal to the Lie algebra element X w . 

Let us now see that it is quite a reasonable thing to do from the point of 
view of the canonical quantum theory. In the case Ai is compact, which is the 
case of interest for us here, a spatial slice £ of M consists of a finite number 
of circles S . The states of the canonical theory are just class functions of the 
holonomies of A, which is the pullback of A on S, around this circles. The 
corresponding field E, which in the quantum theory becomes the operator 
of variational derivative with respect to A. has thus a constant norm along 
each disjoint component of S. It thus makes sense to replace the E field by 
a collection of Lie algebra elements X w constant on each wedge. Moreover, 



as we shall see, the integration over group elements g e in (22) renders all 
X w to have the same norm in each disjoint component of «M, which is in 
agreement with what we get in the canonical approach. 

The group elements g € in (^) are "shared" by two dual faces. Thus, all 
integrals over g e in (|22| ) have as the integrand the product of two matrix 
elements of g e . These can be taken using the formula ( |142D of the Appendix. 
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Figure 6: In two dimensions each face is a union of three wedges, and 
Jit J2, J3 are the currents J w corresponding to these wedges. 

To describe the result of this integration we will assume that A4 has a single 
connected component. Then the generating functional in two dimensions is 
given by the following expression 



Z 2 {J^) = Y J {^ j ) V - E T{P{Ji)P{J2)P{Jz)[ ® (25) 

3 f 




In the integration over the group elements g e in (22) survive only the terms 
in which all spins j a are equal. This spin is denoted by j in (25). The 
symbols V, E stand in ( ^5|) for number of vertices and number of edges in 
the triangulation correspondingly. The product here is taken over all centers 
of faces of A, or, equivalently, over faces /. The graphical notation stands 
for 



Xj(exp Jiexp J 2 exp J 3 ), (26) 



where J\ , J2 , J3 are the three Lie algebra elements (currents) "shared" by 
the face /. Recall that in two dimensions each face is just a union of three 
wedges, and Ji, J2, J3 are the currents J w corresponding to these wedges (see 
Fig. 6). The expression (|25| ) is the final result for the generating functional 
in two dimensions. 

Note that the graph in this graphical representation of the character can 
be obtained as a result of the following simple construction. Let us draw a 
circle S 1 centered at the center of the face (see Fig. 6). The intersection of 
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Figure 7: Tetrahedron from a triangulation of a 3-dimensional spacetime 
manifold. The figure also shows the wedges that lie inside this tetrahedron. 

wedges with the circle gives the circle itself. The intersection of this circle 
with the dual edges gives three points on the circle. Each wedge is labelled by 
a spin j, which is the same for all wedges, and we can assign a label j to the 
circle that intersects the three wedges. Then the face contribution to (p5|), 
that is, the contribution that we graphically represent by (|2.2j), is simply 
the spin network constructed with the labelled circle, with insertion of the 
three group elements (exp Ji exp J 2 exp J 3 ) at the points where the circle is 
intersected by the corresponding wedges. This construction is trivial in two 
dimensions, but turns out to be generalizable to any dimension. 

2.3 3 dimensions 

The three-dimensional case is more interesting. Here the dual faces no longer 
cover the spacetime manifold, and the E field acquires a true distributional 
character. The general construction prescribes to replace the E field, which is 
now a one-form, by a distribution concentrated along dual faces and constant 
along the wedges. It is illustrative to give a coordinate expression for such a 
distributional field. Let us consider an arbitrary wedge w (see Fig. 7). Let 
u = be the equation of a plane containing w. Then the field E that has 
the correct distributional character is given by the expression 

du5(u)X w . (27) 



To calculate the generating functional let us note that each edge of the 
dual triangulation belongs to three dual faces. Thus, in (22), one has to 
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take integrals over products of three matrix elements. The corresponding 
formula is given in the Appendix The result of this integration can be 
described as a generalization of the "circle" construction given at the end of 
the previous subsection. Let us draw a sphere S 2 centered at the center of 
each tetrahedron. The wedges belonging to a particular tetrahedron intersect 
the sphere and draw a graph on its surface. We will refer to this graph as 
Tt, where t refers to a tetrahedron of the triangulation that was used to 
construct Tt- It is not hard to see that Tt is a tetrahedron whose vertices 
come from the intersection of the sphere with the dual edges. The same 
graph Tt can also be obtained by looking at the boundary of tetrahedron t 
(this boundary has the topology of S 2 ), which is triangulated by the faces 
of t, and constructing the graph dual to the triangulation of this S 2 . The 
resulting tetrahedron is the same as that in the previous construction with 
a sphere. 

In (^2|) the sum is taken over spins j a labelling the dual faces. The result 
of the integration over the group elements g e will still be a sum over spins 
j a . Each wedge belongs to some dual face, and, thus, is labelled by spin. 
As we have just seen the edges of Tt are in one-to-one correspondence with 
the wedges. Let us label these edges with the same spins as those of the 
corresponding wedges. We have a current J w associated with each of the 6 
wedges belonging to tetrahedron t. Let us denote this currents by J%, . . . , 
We can now construct a function of the group elements exp J\ , . . . , exp Jq 
that is just a spin network function. 

The generating functional in three dimensions is then given by the fol- 
lowing expression 




where the sum is over all possible coloring of dual faces. This is the final 
expression we are going to use in the following section in order to derive spin 
foam models. 
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2.4 4 dimensions 

The case of four dimensions is analogous to the just analyzed case of three 
dimensions. The only difference is that it is harder to visualize a four- 
dimensional triangulated manifold, and that the final expression for the 
generating functional is more complicated. However, the final result for 
the generating functional follows the same pattern. 

First, let us give a coordinate expression for the distributional field E 
in four dimensions. Let u, v be two functions such that u = 0, v = is 
the equation for a 2-surface in A4 containing one of the faces a of the dual 
complex, and duA dv is positive as defined by the orientation of a. We then 
define E to be equal to 



on each wedge w, and to be equal to zero everywhere else. We repeat this 
construction for all wedges, and add all these distributional forms together 
to get a two-form that is concentrated along the faces of the dual complex. 
Then the integral of such distributional E over a wedge of the triangulation 
is equal simply to X w corresponding to that wedge. 

Let us now describe the result of integration over the group elements 
g e in (^). In four dimensions each dual edge is shared by four dual faces. 
Thus, one has to take the integral of the product of four matrix elements. 
The required formula is given in Appendix ^. The result of this integration 
can again be described using a certain spin network function of the group 
elements corresponding to the currents J w . As in the case of two and three 
dimensions, let us introduce special graphs in the vicinity of the center of 
each 4-simplex h, which we will call T^. We define graph as the intersec- 
tion of the sphere 5 3 surrounding the center of 4-simplex h with the dual 
faces a. The graph lives in S 3 , but it can be projected on a plane. This 
gives us the pentagon graph. Note that at this stage we do not care about 
types of crossings we get (that is, whether this is under- or over-crossing). 
The edges of this graph are in one-to-one correspondence with the wedges w 
belonging to the 4-simplex h. Thus, we can associate to each edge a current 
J w , and label it with the spin j CT labelling the dual face to which the corre- 
sponding wedge belongs. Additionally, let us label each of five vertices of the 
pentagon by a half-integer (spin). Vertices of the pentagon are in one-to-one 
correspondence with tetrahedra of the triangulation. Thus, we shall use the 
notation jt for this spins. 

Inside each four-simplex there are 10 wedges. Let us denote the corre- 
sponding currents by J\ , . . . , Jio . One can then construct the function of the 
group elements exp J\ , . . . , exp J\q . 




(29) 
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The generating functional Z4 is given by the sum over spins of products 
over 4-simplices of the above functions of the currents J w 



Zi(J, A) = £ J] dim,. J] dim Jt J] P{.h) ■ ■ ■ P(J 10 ) 

3<r,jt °" t h jf 

(30) 

Here h labels 4-simplices of the triangulation and j a are spins labelling the 
dual faces. The resolution of each vertex in this pentagon graph is given in 
the Appendix |^. 



3 Theories of interest 

Having calculated the generating functional in various spacetime dimensions 
we are in the position to use our results to study concrete theories. However, 
let us first specify in more details the class of theories that can be treated 
within our approach, and give physically interesting examples of theories 
belonging to this class. In this paper we will restrict our attention to a special 
class of theories. For simplicity we fix the gauge group to be G = SU(2). Let 
A4 be an oriented smooth d-dimensional manifold, and let P be an SU(2)- 
bundle over Ai . The basic fields of a theory of interest are a connection A on 
P, and an Ad-P-valued (d-2)-form, often called B, which we call E because 
of its relation with the SU(2) "electric" field of the canonical formulation of 
the theory. The action of the theory is of the form 

S[A,E]= [ [Tr(E A F) + $(£)] , (31) 
JM 

where F is the curvature of A, 'Tr' is the trace in the fundamental rep- 
resentation of SU(2), and $ is a gauge invariant polynomial function that 
depends on E but not on its derivatives. Also, &{E) should be a d-form 
so that it can be integrated. It can also depend on other dynamical fields, 
for example Lagrange multipliers as in the case of gravity, or on additional 
non-dynamical fields (background structure), as in the case of Yang-Mills 
theories. Let us now give examples of theories belonging to this class. 
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3.1 BF theories 



The action of BF theory in any dimension is given by 




JM 



that is, this is the action of the type (|3l| ) with zero "interaction" term 
The equations of motion that follow from this action state that 



d A E = 



F = 



(33) 



where d^E is the covariant derivative of E. Thus, the flat connections play 
the dominant role in BF theory. 



Yang-Mills action can be written only if one introduces a background metric 
on Ai. However, in two spacetime dimensions, Yang-Mills action turns out 
to depend only on a measure (area two- form) defined by the metric on Ai. 
Thus, we will keep track only of this dependence of the action on a measure 
djj, on Ai. Yang-Mills theory in 2D is described by the following BF-like 
action 



Here E is a Lie-algebra valued zero-form, F is the curvature of the con- 
nection A. Solving equations of motion for E that follow from this action, 
and substituting the solution into the action, one can check that the action 
reduces to the standard Yang-Mills action 



where g a b is a metric (of Euclidean signature) on Ai such that d 2 Xyfg = dfj,, 
where y/g is the square root of the determinant of the metric. 

Using the techniques developed in this paper, one can calculate the 
vacuum-vacuum transition amplitude of this theory, that is the path in- 
tegral of exp {i x action). However, an interesting feature of this theory is 
that the same techniques can be used to calculate the partition function, 
that is, the path integral 



3.2 2D Yang-Mills 




(34) 




(35) 




(36) 
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The problem of calculation of the partition function can be reduced to the 
problem of calculation of the vacuum- vacuum transition amplitude of a some- 
what different theory, that is of the path integral of exp (i x action), with 
the action given by 

/ Tr (EF) — i%- I dfiTr(E 2 ). (37) 
JM 2 Jm 

Indeed, integrating over E in the path integral of exp (i x action), one gets 
(|36[). As one can see, the two actions (|3~4T) , (|37|) differ by the factor of i in 
front of the second term. Using the two actions (|34"1) , (|37|) we will calculate 
below both the vacuum- vacuum transition amplitude and partition function 
of the theory. 



3.3 3D BF theory with the 'cosmological term' 

In three dimensions the E field of BF theory is a one-form, and one can add 
a term cubic in E to the BF action to obtain 

- / Tr (e A F + —E AEAe], (38) 
J M V 12 / 

where j\A is assumed to be a three-dimensional orientable manifold. The 
action is a functional of an SU(2) connection A, whose curvature form is 
denoted by F, and a 1-form E, which takes values in the Lie algebra of 
SU(2). Thus, the action (^) is that of BF theory in 3d, with E field playing 
the role of B, and with an additional "cosmological term" added to the usual 
BF action. This theory is related to gravity in 3D as follows. Having the 
one-form E, one can construct from it a real metric of Euclidean signature 

9ab = -~Tr(E a E b ). (39) 

Here o, b stand for spacetime indices: a, b, ... = 1, 2, 3. We take the E field to 
be anti-hermitian, which explains the minus sign in (|39|). Thus, the E field in 
( |38|) has the interpretation of the triad field. One of the equations of motion 
that follows from ( |3"8| ) states that A is the spin connection compatible with 
the triad E. Taking the E field to be non-degenerate and "right-handed", 
i.e., giving a positive-definite volume form 

±s abc Tr(E a E b E c ), (40) 

and substituting into (|38|) the spin connection instead of A, one gets the 
Euclidean Einstein-Hilbert action 

\ f d 3 x^(R-2A). (41) 
2 Jm 
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We use units in which SttG = 1. The minus in front of (^) is needed to yield 
precisely the Einstein-Hilbert action after the elimination of A. Thus, A in 
(|38|) is the cosmological constant. An important difference of this theory 
and Einstein's theory is the fact BF theory, unlike gravity, is defined for 
degenerate metrics. 



3.4 4D BF theory with the 'cosmological term' 

In four spacetime dimensions the E field of BF theory is a two-form. Thus, 
one can construct the following action functional 

/ Ti(E A F) — —Tr(E A E), (42) 

where A is a real parameter, which we will refer to as the "cosmological con- 
stant" . If one, as in the following subsection, adds to this action an additional 
constraint that E is simple, that is, given by a product of two one-forms, 
then this theory is equivalent to Einstein's theory and A is proportional to 
the "physical" cosmological constant. 



3.5 4D Euclidean self-dual gravity 

The action for Euclidean general relativity in the self-dual first order for- 
malism is given by [24] 

J Tr(E A F) — ^ (e? A E^ - h lj E k A E^j , (43) 

where, to write E l , we have introduced a basis in the Lie algebra of SU(2). 
Here tp 11 is a symmetric matrix of Lagrange multipliers. The variation of 
the action with respect to ip yields equations 

E l AE j = 5 tj -E k AE k . (44) 
3 

In the case E is non degenerate, i.e., the right hand side of ([i4|) is non 
zero, these equations are satisfied if and only if E is the self-dual part of a 
decomposable 2-form, i.e., if and only if there exists a tetrad field e 1 , I = 
0,1,2,3 such that E % = ±(e° A e l + \eij^ A e k ). In this case the action 
reduces to the self-dual Hilbert-Palatini action for 4-dimensional gravity and 
ipij correspond to the components of the self-dual part of the Weyl curvature 
tensor. Thus, (^) is the BF action in four dimensions, with an additional 
"simplicity" constraint added to it. 
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4 Corresponding spin foam models 

In this section we present state sum models for the theories we just described. 



4.1 Ponzano-Regge-Ooguri models 

These models give an expression for the vacuum-vacuum transition ampli- 
tude of BF theory, i.e., the path integral 

Z BF {M) = J VAVEexp (i J Tr(E A F)^j , (45) 

as a sum over colored triangulations. More precisely, let us fix a triangulation 
A of M. In three dimensions, let us label the edges of A by irreducible 
representations of SU(2), that is, by spins j. Thus, to each edge e we assign 
a label j e . One can then construct the following sum over labellings 

pR(A)=^n dim >iK 6 ^- ( 46 ) 

Here dim(j) = 2j + 1 is the dimension of the representation j, the second 
product is taken over tetrahedra t of A, and (6j) is the (normalized) classical 
(6j)-symbol (see the Appendix C for a definition) constructed from the six 



spins labelling the edges of t. Note that the sum over spins in (46) diverges. 
To make sense of it one must introduce a regularization. A possible regular- 
ization is given by the Turaev-Viro model, which we discuss below. After the 
introduction of this "regularization" , this sum can be shown to be triangu- 
lation independent. Thus, (flG|) gives an invariant of A4: PR(A) = PR(A^). 

In four dimensions one can construct a similar sum by labelling the faces 
of A. Let us denote the spin labelling a face / by jf. Let us consider the 
following quantity: 

o(A) = e n dim ^ n dim Jt n( i5 i)^ ^) 

jf,jt f t h 

After a certain "regularization" (given by the Crane- Yetter model below) 
that gives sense to the infinite sums in 

(H), this 

can be shown to be trian- 
gulation independent: 0(A) = 0(M). In ([17]) one has introduced an addi- 
tional label jt for each tetrahedron of A. The spin jt labels an intertwiner 
one has to assign to each tetrahedron (see the subsection on Crane- Yetter 
model for details). The last product is taken over the 4-simplices h of A. 
Then (15 j) is the (normalized) (15j)-symbol constructed from ten spins jf la- 
belling the faces of h and five spins jt labelling the five tetrahedra composing 
h. See the Appendix |C| for a definition of the normalized (15j)-symbol. 



L. FREIDEL, K. KRASNOV 



1211 



As was realized by Ooguri Q, these models give a "discrete" realization 
of the path integral because they can be obtained from the requirement that 
the connection on M is flat. Indeed, formally the path integral (^5|) is equal 
to 

r VA5(F). (48) 



The state sums ( [46|) , ( f47|) are in certain precise sense realizations of the inte- 
gral over connections with the integrand being the delta-function at F = 0. 



4.2 Migdal-Witten model 

Migdal || has studied the Yang-Mills theory on a lattice, and, in particular, 
has proposed a lattice model for the Yang-Mills theory in two dimensions. 
This model was later studied by Witten Q in the connection with topological 
field theories in two dimensions. 

As we discussed above, the problem of calculation of the Yang-Mills 
partition function can be reduced to the problem of calculation of the path 
integral 

Z YM {pe 2 ,M) = J VAVEexp (^i J^Tr(EF) + & — dpTr(E 2 )j , (49) 

where p is the total area of A4, and other notations are explained in sub- 
section (|3.2j ). Migdal ||| has proposed the following lattice version of this 
path integral. For simplicity, we will formulate Migdal model on a trian- 
gulated manifold. Let us triangulate our two-dimensional manifold M, and 
introduce the dual triangulation (see Fig. 5). Let us label the edges e of 
the dual complex by group elements g e , and the dual faces a by irreducible 
representations of SU(2), i.e., spins j a . The "discrete" version of (|49|) is then 
given by (see ||): 

YM(pe 2 ,.M) = [] / ^En dim >^(^i---^Jexp(-eVc(i CT )/2) , 

(50) 

where the multiple integral is performed over all group elements g e and dg 
is the normalized Haar measure on SU(2), p a is the area of the dual face 
a, as defined by the measure dp, c(j) is the value of the quadratic Casimir 
operator in the representation j 

c(j) = 2j(j + 1), (51) 

and g ei ■ ■ ■ g en is the product of group elements around the dual face. Af- 
ter integration over the group elements g e the partition function takes the 
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following simple form: 

3 

where k(M) is the Euler characteristics of M. Thus, as we indicated in 
the argument of YM, the partition function depends only on the topological 
properties of M. Also, the dependence on measure d\i enters only through 
the dependence on the total area p of M.. 



4.3 Turaev-Viro model 



The Turaev-Viro model gives a way to calculate the transition amplitude of 
the theory defined by (j38[), i.e., the path integral 



Zf F (A,M) = JvEVAexp(^-i J Tr (e A F + A E A E\\ (53) 

In this paper we consider only vacuum-vacuum amplitudes. Although the 
Turaev-Viro model can be used to calculate more general amplitudes between 
non-trivial initial and final states, we will not use this aspect of the model 
here. We consider the version of the model formulated on a triangulated 
manifold. Thus, let us fix a triangulation A of Ai. Let us label the edges, 
for which we will employ the notation e, by irreducible representations of 
the quantum group (SU(2)) 9 , where q is a root of unity 



2-1 



q = e~ = e in . (54) 

Later we will relate the parameter h with the cosmological constant A. The 
irreducible representations of (SU(2)) 9 are labelled by half- integers (spins) 
j satisfying j < (k — 2)/2. Thus, we associate a spin j e to each edge e. 
The vacuum-vacuum transition amplitude of the theory is then given by the 
following expression (see, for example, p5f| ): 

TV(g, A) = r, 2V £ ndim,(j e ) \\^j) q , (55) 

je e t 

where r/ and the so-called quantum dimension dim g (j) are defined in the 



Appendix A by ( |129| ) and ( 130 ) correspondingly, and V is the number of 
vertices in A. The last product in ( |55| ) is taken over tetrahedra t of A, 
and (6j)q is the (normalized) quantum (6j)-symbol constructed from the 6 
spins labelling the edges of t (see Appendix C). It turns out that (|55| ) is 
independent of the triangulation A and gives a topological invariant of A4: 
TV(g,A) = TV(g,M). 
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The construction that interprets the Turaev-Viro invariant ( |55|) as the 
vacuum-vacuum transition amplitude of the theory defined by (pq) is as 
follows. It has been proved (see e.g. pER ) that ( |55| ) is equal to the squared 
absolute value of the Chern-Simons amplitude 

TV(q,M) = \CS(k,M)\ 2 , (56) 



with the level of Chern-Simons theory being equal to k from (54). It is 



known, however, that the action (38) can be written as a difference of two 
copies of Chern-Simons action 

S CS (A) = ^- [ Tr ( A AdA + -A A A A a] . (57) 
47r Jm \ 3 J 



Indeed, note that 

k\ r f A 2 \ 

S(A + XE)-S(A-XE) = — Tr [E A F + —E A E A E , (58) 

n Jm \ 3 J 

where A is a real parameter. Thus, (58) is equal to ( |3lf ) if 

A = -VA/2, k = -=, or n = VA. (59) 
VA 

This relates the deformation parameter q of the Turaev-Viro model with 
the cosmological constant A, and proves that the Turaev-Viro amplitude 
is proportional to the vacuum-vacuum transition amplitude of the theory 
defined by @ 

TY(q,M) oc Z 3 BF (A,7W). (60) 

To compare in Sec. | the spin foam model obtained via our techniques 
with the Turaev-Viro model, we will need the first-order term in the de- 
composition of (^) in the power series in A. This is proportional to the 
expectation value of the spacetime volume in BF theory at A = 0. Indeed, 
the expectation value of the volume is given simply by the derivative of the 
amplitude (0) with respect to (— iA) evaluated at A = 0: 

( v„i, = fv TJZ°] ( " )e ' s = ■ ( fl '"!f (A) ) , (si) 

/ A=0 



/ VEVA e iS \ <9A 

Vol(.M) = / ±-e abc Tr(E a E b E c ). 
Jm 1^ 



Thus, the expectation value of the volume of M in BF theory is proportional 
to the first order term in A in the decomposition of (55). 
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Let us find this expectation value. An important subtlety arises here. It 
is not hard to show that ( f35l) has the following asymptotic expansion in h 



l^j PR(A) (l - ft 2 i(Vol>) , (62) 

where PR is the amplitude of the Ponzano-Regge model (|4q) V is the number 
of vertices in A, and i(Vol) is a real quantity independent of h. Thus, appar- 
ently there is no term proportional to H 2 = A in this expansion. However, 



as it is explained, for example, in pq| , the state sum invariant (55) does not 
exactly give the transition amplitude of BF theory. Instead, the Turaev-Viro 
amplitude is only proportional to the BF amplitude, and the proportionality 
coefficient depends on h. This can be understood as follows. The integration 
over (A + XE), (A - \E), which is carried out to obtain \CS(k,M)\ 2 in (||) 
and thus the Turaev-Viro amplitude, is different from the integration over 
A, E one has to perform to obtain (|53]). The difference in the integration 
measures is a power of U. Thus, the amplitude (|53|) and the squared abso- 
lute value of the amplitude of the Chern-Simons theory are proportional to 
each other with the coefficient of proportionality being a power of fi. In the 
discretized version of the theory, given by the Turaev-Viro model, this power 
of h is replaced by h 3V . Thus, the Turaev-Viro amplitude, in the limit of 
small cosmological constant, differs from the BF amplitude by a power of 
h 3V . 

This remark being made, we can write an expression for the expectation 
value of the volume in BF theory: 

7V1\ 3 ( TV ( A > A ) \ 

' ^iVol(A,j)fn dim (j'e)Il( 6 i))' (63) 



PR(A) 



Je 



t 



where the function Vol(A, j) of the triangulation A and the labels j = {j e } 
is given by 




dA V(&74*)77 A= o e V OA 
dln((6j) q )_ 

A=0 



dA 



A=0 
(64) 



Here v stands for vertices of A, e stands for edges and t stands for tetrahedra. 
We intentionally wrote the expectation value of the volume in the form (|63]) 



L. FREIDEL, K. KRASNOV 



1215 



to introduce the volume Vol(A,j) of a labelled triangulation, which will be 
of interest to us in what follows. Indeed, ( |63| ) has the form 

E je iVol(A,j)Amplitude(A,j) 

E je Amplitude(A,j) ' { > 

where 

Amplitude(A, j) = [] dim(j e ) JJ(6j) (66) 

e t 

is the amplitude of Ponzano-Regge model. This shows that Vol(A, j) indeed 
has the interpretation of the volume of a labelled triangulation. Note that 
the volume turns out to be purely imaginary. This has to do with the fact 
that in (^3|) one sums over configurations of E of both positive and negative 
volume. A more detailed explanation of this is given in p6| . 

The volume (^) has three types of contributions: (i) from vertices; (ii) 
from edges; (hi) from tetrahedra. It is not hard to calculate the first two 
types of them. One finds that each vertex contributes exactly 1/12, and 
each edge contributes j e {je + l)/6, where j e is the spin that labels the edge 
e. It is much more complicated to find the tetrahedron contribution to the 
volume, that is, the derivative of ln((6j) 9 ) with respect to A. The result is 
described in the Appendix D. Using the notations of the Appendix D, the 
final result for the expectation value of the spacetime volume in BF theory 
can be written as follows: 




Here ^ stands for the normalized classical (6j)-symbol. 



It is interesting to note that not only tetrahedra t of A contribute to 
the volume, but also the edges e and vertices v. The contribution from the 
vertices is somewhat trivial - it is constant for each vertex. Nevertheless, 
when thinking about the triangulated manifold A4 one is forced to assign 
the spacetime volume to every vertex. The contribution from edges depends 
on the spins labelling the edges. Again, this implies that each edge of the 
triangulation A carries an intrinsic volume that depends on its spin. The 
contribution from tetrahedra is more complicated. It is given by a func- 
tion that depends on the spins labelling the edges of each tetrahedron. It is 
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interesting that this picture of the spacetime volume being split into contri- 
butions from vertices, edges and tetrahedra can be understood in terms of 
Heegard splitting of A4 . Recall that Heegard splitting of a three-dimensional 
manifold A4 decomposes A4 into three dimensional manifolds with bound- 
aries. Then the original manifold can be obtained by gluing these manifolds 
along the boundaries. For the case of a triangulated manifold A4 , as we have 
now, the Heegard splitting proceeds as follows. First, one constructs balls 
centered at the vertices of A. Then one connects these balls with cylinders, 
whose axes of cylindrical symmetry coincide with the edges of A. Removing 
from A4 the obtained balls and cylinders, one obtains a three-dimensional 
manifold with a complicated boundary. One has to further cut this mani- 
fold along the faces of A. One obtains three types of "building blocks" that 
are needed to reconstruct the original manifold: (i) balls; (ii) cylinders; (hi) 
spheres with four discs removed. Each of this manifolds carries a part of the 
original volume of A4. Our result ( |67| ) provides one with exactly the same 
picture: the volume of A4 is concentrated in vertices (balls of the Heegard 
splitting), edges (cylinders), and tetrahedra (4- holed spheres). 



4.4 Crane- Yetter model 



Crane, Kauffman and Yetter [10] studied a state sum model that is very 
similar to the Ooguri model |J, however, instead of the gauge group SU(2) 
the quantum group (SU(2)) 9 is used. The state sum invariant they proposed 
is given by (see, for example, [pij|]): 

CY(q, A) = r]- 2V+2E E E II dim ^/) II dim <^) IK 15 iV ( 68 ) 
jj it f t h 

Here dim g (j) is the quantum dimension (|130|) , 77 is defined by (|12S| ), (15j) g 
stands for a (normalized) quantum (15j)-symbol (see the Appendix |C|) as- 
sociate with any 4-simplex. This quantum (15j)-symbol is the Reshetikhin 



Turaev evaluation of the graph described in subsection (2.4). The symbols 
V and E stand for the number of vertices and edges in A correspondingly. 
This state sum model is independent of a triangulation of the 4-manifold 
used to compute (^8|). Thus, it gives an invariant associated with the mani- 
fold M: CY(q,A) =CY(q,M). 

The quantity CY(q,A4) is conjectured to give the transition amplitude 
for BE theory with cosmological constant, the cosmological constant A being 
related to q by 

q = exp iA. (69) 

Thus, unlike in the case of 3-dimensions, we have now h = A. This relation 
can be established as follows. It has been shown (see e.g. pl| ) that if A4 is a 



L. FREIDEL, K. KRASNOV 



1217 



manifold with boundary then CY(g, M) is proportional to the Chern-Simons 
transition amplitude CS(A:, dM) introduced in the previous subsection: 

CY(q,M) oc CS{k,dM). (70) 

The deformation parameter q of the Crane- Yetter model is related to the level 
of Chern-Simons theory as in fl54|). Consider now the transition amplitude 
for BF theory with cosmological term: 

Zf F (A,M) = J VAVEexp (i J Tr(E A F) - ^Tr(£ A E)^j . (71) 

Let us integrate over the field E in this path integral. One gets 

Z 4 BF (A,.M)oc |pAe X p^^Tr(FAF)). (72) 

Now, using the fact that 



Tr(F A F) = d(A A dA + -A A A A A) (73) 

3 



we get 



Zf F (A,M) oc J VAexp (±-S cs (A, dM)\ . (74) 

This, together with (|70|), then implies that 

CY(q,M) oc Zf F (A,M), (75) 

where q is related to A as in (pG|), 

To compare the spin foam model we obtain below with the Crane- Yetter 
model we will need to know the first order term in the decomposition of 
Zf F (A, A4) in the power series in A. To find it, we have to take into account 
the fact that, similarly to the case of 3D BF theory, the transition amplitude 
is only proportional to the Crane- Yetter amplitude (|68|), with the propor- 
tionality coefficient being a function of Ti. This proportionality coefficient is 
given by a power of %. As one can check, the "zeroth" order term in h of 
is equal to 

(h 



0(A), (76) 

where O(A) is the Ooguri state sum invariant (|47|). Thus, the first order 
term in A of Zf F (A,Ai) is given by 
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However, unlike the case of 3D BF theory with cosmological term, in this 
case the only contribution to this expression comes from the quantum (15j)- 
symbol. All other "building blocks" that one encounters in (65) contain only 
the h 2 = A 2 terms, and, thus, do not contribute to (|77|). The first derivative 
of the quantum (15j)-symbol with respect to h = A can be found by methods 
analogous to the ones used in the Appendix D to calculate the derivative of 
(6j)-symbol. However, in the case of (15j)-symbol, the calculation is much 
simpler due to the fact that only crossings contribute to the first order in Ti. 
The (15j)-symbol used in (|68| ) contains only one crossing (see, for example, 
p5|]). Using the expression ( |152| ) for the R-matrix as a formal power series 
in h given in the Appendix D, one can easily check that the first border 
term of the quantum (15j)-symbol is given by: 




where the dots stands for terms containing graspings between two edges 
sharing one vertex and correspond to different possible choices of the framing 
of the 15j-symbol. We will not keep track of these terms at this stage. Their 
relevance will be discussed below. This is, however, not quite the expression 
we want because it is not symmetrical. We will symmetrize it by putting 
the grasping at all pairs of lines of the graph in (78) that do not share a 
vertex. However, if one does that, than the quantity one obtains for each 
given 4-simplex h is not equal to ( f78|) . Only if one sums over all 4-simplices 
of the triangulation one obtains the quantity that is equal to the sum of ( |Tq) 
over h. This can be proved by using certain simple relations that hold for 
graspings. One such relation is the analogs of "closure" relations ( |160| ). It 
graphical representation is given by: 




y = (79) 



Such "closure" relation holds for every vertex of the graph in (ff^). Another 
relation that one has to use involves two different pentagon graphs. As one 
can convince oneself, the vertices of the pentagon graph in (|78]) are in one- 
to-one correspondence with the tetrahedra of A. Thus, every tetrahedron is 
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shared exactly by two 4-simplices, and the following relations holds: 







= 



(80) 



This allows one to cancel certain types of graspings from one 4-simplex 
with similar types of graspings from the neighboring 4-simplex. Certain 
cancellations occur for each tetrahedron of the triangulation. These relations 
allow one to obtain the following symmetric expression: 



/ 



1 x - 1 

2 V 30 

h 



^sign(e,e') 



e,e' 



V 



e' 




(81) 



where the sum is taken over different pairs e, e' of edges that do not share 
a vertex. There are 30 terms in this sum, which explains the factor of 1/30 
in (|8l|). The quantity sign(e, e') in (|8l|) is plus or minus one depending on 
the orientation of the two edges e,e'. A consistent choice of orientations 
comes from the geometrical 4-simplex discussed in the Appendix E. With 
this choice, sign(e, e') = sign(/, /'), where we use the fact that every edge 
of the pentagon graph in (|8l| ) is in one-to-one correspondence with a face of 



the corresponding 4-simplex, and sign(/, /') is defined by the equation ( 161 ). 
The dots in (81) correspond to some symmetric expression that contains only 
terms with graspings between two edges sharing one vertex. 



4.5 Reisenberger model 

Reisenberger || has proposed a state sum model corresponding to the self- 
dual Plebanski action ([1^) . The main idea of his construction is to modify 
the SU(2) Crane- Yetter model in such a way that the constraints ( f44| ) are 
implemented. The model can be described as follows. As in the previous 
subsection, let M. be a triangulated 4-manifold. Let us associate a group 
element g w to each wedge w of the dual triangulation (see Section |2[) . One 
can think of g w as the holonomy of connection around the boundary of the 
wedge w with the basepoint being the center of the 4-simplex h to which the 
wedge belongs. Let D l w be a differential su(2) operator acting on functions 
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of g w as the sum of left and right invariant vector fields. Let us denote by 
O lJ the following operator acting in (SU(2))® 10 : 

W = £ e(w, w')DlDl - \^ £ e(w, w')D k w D w , k , (82) 

w,w' w,w' 

where e(w,w') is the sign of the four volume span by the wedges w, w' 
(e(w,w') = if the two wedges don't span a four-volume). We choose a 
coloring of the wedges by spins j w , denote by Rj(g) the representation of 
g as an endomorphism of Vj, and define for each 4-simplex h the following 
function on (SU(2))® 10 : 

<t>h,jig w ) = tr® w&hV]w ^-j=e~^ mQii ®w&hv 3w (2i w + l)Rj w (9w)j ■ (83) 
The state sum model proposed by Reisenberger is given by 

So ? / II d 9™ II <f>hj(gv, ) (84) 



J w 



The limit is a possible way of selecting only the states belonging to the kernel 
of Wi. 



5 Applications of the generating functional 

In this section we use the generating functional computed in Sec.|2| to con- 
struct systematically state sum models of the theories described in Sec.|3[ 
We then compare the results of this systematic approach with the known 
results described in Sec.|I[ 



5.1 BF theory: any dimension 

The transition amplitude of the theory is simply the value of the generating 
function Z at J = 0. As one can easily see, this gives exactly the Ponzano- 
Regge model (46) in three dimensions, and the Ooguri model (E7|) in four 



dimensions. The corresponding two dimensional model is somewhat trivial: 
the transition amplitude is given by (p5|) at J = 0, which gives 

J2(d^j) xlM \ (85) 

3 

where (M) is the Euler characteristics of M. One can recognize in this 
expression the volume of the space of flat connections on A4 modulo gauge 
transformations expressed in term of Riemann zeta function (see, for exam- 
pie, p|). 
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5.2 2D Yang-Mills 

Our general strategy is to calculate the 'interaction' term &{E) of the action 
( |3l|) on the configuration of E that is distributional along the wedges of the 
dual complex 

E = Y,E W , (86) 

w 

express the result as a polynomial function &(X) of the variables X w , in- 
troduced in the previous section, and then look for the vacuum-vacuum 
transition amplitude of the theory as given by 

(e&l-VW Z[J,A]) (87) 

In the case of 2d Yang-Mills theory, the 'interaction' term i<3? is given by 
(see©) 

|/ ^Tr(£ 2 ), 
2 Jm 

Calculating this on the configuration of E given by (^6|), with each E w being 
constant along the wedge w, we get 



e 2 



5>„,Tr(X 2 ) = -e^P^KK, 



where p w is the area of the wedge w, as measured with respect to dfi, and 
we have introduced the SO(3) indices ( |128| ). The partition function of the 
theory is then given by 



cxp 



J=0 



To compare this with the partition function of the Migdal-Witten model, 
it is more convenient to take the generating functional Z[J] in its general 
form (|2|). Then, using 



e Q G4 ,4) . P(J)R {j) ( e J )\j =0 = e a ^ +1 ^ 2 = e «Mi)/2+i/4) (91) 

where R(j\(e ) is the group element e J taken in the j representation, and 
c(j) is the quadratic casimir of this representation. We get for the partition 
function 

II / d ^ EIlH^fei • • • 9eJ exp j -e 2 p w {c{j a )/2 + 1/4) J . 

(92) 
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This can be seen to be equal to 

YM{pe 2 ,M)e- e2p/A (93) 

by noting that 

^2 Pw = Pa, and Y^Pa = P- 

wda a 

Thus, our approach gives the correct expression ( |50j ) for the partition func- 
tion of the Yang- Mills theory, apart from the factor exp(— e 2 p/4). However, 
the later is what is called a "standard renormalization" of the partition 
function. That is, the result for the partition function of 2D Yang-Mills may 
depend on a regularization procedure that was used to calculate it, the two 
different schemes giving results that may differ, in particular, by factors of 
exp(— ae 2 p), where a is some coefficient. See, for instance, |1 for a more 
detailed discussion of the standard renormalizations. Thus, the partition 
function obtained by our method differs from that of the Migdal-Witten 
model just by a standard renormalization factor. 



5.3 3D BF theory with cosmological term 



As one can see from (38), the 'interaction' term i<&(E) for this theory is 
given by 

i$(E) = -— [ Tt{EAEAE). (94) 
12 Jm 

As in the case of Yang-Mills theory in two dimensions, one has to find a 



polynomial function &(X) of variables X w by evaluating (94) on the distri- 
butional field E given by (||) . Unfortunately, in the case of three dimensions 
the result is not well-defined as in the case of 2D Yang-Mills theory. Indeed, 
the E field is concentrated along the wedges (see Fig. 7) , and non-trivial con- 
tributions to the integral ( |94j ) come from the points where wedges intersect. 
Thus, the contributions to (|94] ) come from the integrals 

/ Tr(E Wl AE W2 AE W3 ), (95) 
Jm 

where 101,102,103 are three wedges that intersect. The wedges can intersect 
at points, as, for example, they do at the center of each tetrahedron, or they 
can have more general intersections, as, for example, an intersection of three 
wedges sharing a dual edge. Or, instead, one can have in (|9q) w\ = W2- 
Whatever the type of an intersection is, the result of (|95|) is ill-defined be- 
cause of the distributional nature of E field. In this paper we show how 
the ambiguity in ( p5|) can be resolved for the simplest, and most important 
type of intersections: intersections of three different wedges at the centers of 
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tetrahedra. As we shall see, for intersections of this type there is a simple 
way to resolve the ambiguity in (|9^) using geometrical considerations. As 
we show, this type of intersection is the most important one, for it is respon- 
sible for, in certain sense, the most intricate contribution to the transition 
amplitude of the theory. Thus, in this subsection, we restrict our consider- 
ation only to this special type of intersection, calculate the corresponding 
interaction term <&(X), find the corresponding spin foam model, and com- 
pare it with the Turaev-Viro model. In the first order in the decomposition 
of Turaev-Viro amplitude in power of A, we will be able to identify a term 
analogous to the term we obtain in our model and compare them. We will 
find a very delicate matching between the two, including the numerical co- 
efficients. The importance of other types of intersections, not treated here, 
will be emphasized later. 

Thus, the intersections we consider are the ones for which three different 
wedges intersect at the center of a tetrahedron. Given a tetrahedron t, 
there are 20 (without counting the permutations) different triples of wedges 
w i 7^ W2 7^ 1^3- Four of these triples do not span a three-volume; thus, 
it is natural at first to take the integral (95) for such triples to be zero. 



These are exactly the triples of wedges that share a line - a part of the dual 
edge. Thus, there are only 16 different triples of wedges (without counting 
the permutations) that contribute to ( |95"| ) for a given tetrahedron. Recall 
now that each E w is given by (^7|). Thus, for any given triple of wedges 
w\ / u>2 7^ W3, the integral (|9~5|) is proportional to 

iATv{X Wl X W2 X W3 ), (96) 

but the proportionality coefficients are not fixed, because of the indetermi- 
nacy in the value of the integral 

5(ui)5(ii2)5(u3)dui A du2 A du^. (97) 

Of course if the wedges were infinite planes without boundaries this integral 
is a well-defined quantity: it is just the intersection number and as such 
should be ±1 depending on the orientations of the wedges. Because of the 
existence of boundary we expect this integral to be a number smaller than 
1. We fix the proportionality coefficient by requiring that 

^jTr(EAEAE), (98) 

where the integral is taken over a tetrahedron t, is equal to the geometrical 
volume of t. For any triple of wedges wi / »2 / w%, the integral (|97|) is 
equal to a -sign(u)i, u>2, ^3), where a is a coefficient that is independent on a 
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triple, and sign(u;i, W2, W3) is plus or minus one depending on the orientation 
of the form du\ A du2 A du^. Thus, @) is equal to 

— a- siga(w 1 ,w 2 ,w 3 )Tr(X Wl X W2 X W3 ). (99) 

Here the sum is taken over triples of different wedges in t, taking into account 
all different permutations of w±, w 2 , 1^3. Taking into account the fact that 
the number of permutations of wi, W2, 1^3 is 6, we get for (|98[): 

l - J2 aTr(X Wl X W2 X W3 ), (100) 

Wi<W2<W3 

where the notation w\ < w 2 < W3 means that the sum is taken over 16 
different triples w\ ^ W2 7^ u>3 such that signal, W2, W3) = 1. 

To relate ( |100| ) to the volume of tetrahedron t, we recall the geometrical 
interpretation of variables X w . They were introduced in the previous section 
as the variables that carry information about the length of edges of the 
triangulation. At this stage it is more convenient to introduce SO(3) indices 
( |128| ), Thus, each X w is characterized by X l w , i = 1, 2, 3. Recall that wedges 
w are in one-to-one correspondence with edges e of the triangulation. Thus, 
let us view each X l w as the vector representing the corresponding edge (each 
edge can be viewed as a vector pointing from one vertex of A to another), and 
the norm squared X l w X l w of X l w as the length squared of this vector. Indeed, 
it is not hard to check that the interpretation of X' l w X' l w as the length of 
the corresponding edge is consistent with the other known facts. Let us 
consider the operator corresponding to X l w X l w . According to our general 
prescription, this quantity is represented by the operator (5/5iJ w ) 2 . We 
have already dealt with this operator in the previous subsection, see (|9l|). 
Its eigenvalue is given just by the half of the Casimir (plus 1/4). Thus, in 
the sense of eigenvalues, we can write 

Xl J Xl J = (j + l/2) 2 , (101) 

where j is the spin from (^) labelling the dual face that contains w. Thus, 
( |101| ) tells us that, in the limit of large spins j, the norm of X w grows as j. 
This is to be compared with the length spectrum of the canonical quantum 
theory 

(length) = yjj(j + 1). (102) 

This also grows as j for large spins. The expression ( |102j ) can be easily 
derived in the context of canonical (loop) quantum gravity in three dimen- 
sions (note that we use units in which 8irG = 1). Another motivation for 
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interpreting the spin j as the length of an edge (for large j) is that it is ex- 
actly this interpretation that must be used to reproduce correctly the Regge 
calculus version of Einstein-Hilbert action in the Ponzano-Regge model of 
quantum gravity ]7j|. Thus, we learn that the interpretation of the norm of 
each vector X w must be that of the length of the corresponding edge of A. 
Having this fixed we can relate ( JlOCj ) to the volume of tetrahedron t. We 
have 

Tr(X wl X W2 X m ) = 2eijkXl J1 Xi !2 X* 3 , (103) 



where we have introduced the SO(3) indices, see ( [128D . Each X^ has the 
interpretation of the vector corresponding to one of the edges of t. Thus, 



(|103| ) is equal to 12V, where V is the volume of t. In (100) we have 16 such 
terms. Thus, it is equal to 

-a ■ 16 • 12V. 

The requirement that (pq) is equal to the volume of t fixes the parameter a 
to be 1/(16 • 6). Thus, finally, we obtain the interaction term i&(X) to be 

~ "e" ^6 ^2 e vkX l Wl X J W2 X* 3 , (104) 

where the sum is taken over 16 terms. To obtain the transition amplitude 
of the theory we have to replace each X l w by the operator S/SiJ^ and act by 
the exponential of ( 104D on the generating functional. Thus, the first order 



term in A in the decomposition of the transition amplitude is given by: 

This is to be compared with (^). We will now show that the most com- 
plicated term in that expression - the term that involves trivalent graspings 



exactly matches our result (105), including the numerical coefficient and 



the sign. To see this we just have to relate the trivalent grasping in ( 



to the cubic operator in (105). As explained in the Appendix D, a single 
grasping in ( |6"7| ) acts by inserting (a 1 j\[2), where a % are the Pauli matrices. 
The operator (5/5iJ l ), when applied only one time, acts by inserting just 
a 1 . Also, taking into account the definition of the trivalent grasping in (|67|), 
one can show that 



4e 



tjk 



_6 6 5_ 

SiJi SiJL 5U% 



(106) 



e 
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So our approach does account for the most intricate part of the 3D transi- 
tion amplitude in the first A-order, reproducing correctly even the numerical 
coefficients. Note, however, that we do not get all of the terms appearing in 
(|67|). The discrepancy we find between the model we obtain and the usual 
state sum model can be understood by comparing the two different expres- 
sions for the spacetime volume: the one obtained within our approach and 
the one obtained within the Turaev-Viro model, see (|67|). The Turaev-Viro 
model tells us that we must associate the spacetime volume not only to 
tetrahedra of the triangulation, but also to edges and vertices. Within our 
approach, only the tetrahedron part of the spacetime volume is accounted 
for. The reason for this was that, when evaluating the cosmological term of 
the action on the distributional E field, we took only the terms which came 
from intersections of different wedges at the centers of tetrahedra. As we dis- 
cussed, even these terms are ill-defined because of the distributional nature 
of the E field. However, the corresponding ambiguities can be successfully 
resolved for these types of intersections by using geometrical considerations, 
and, after the ambiguities are eliminated, the result exactly matches the 
analogous terms of the usual state sum models, including the matching of 
the numerical factors. The geometrical considerations we used were exactly 
the ones that relate the terms we considered to the volume of a geometri- 
cal 3-simplex. Thus, it is not very surprising that only the "part" of the 
Turaev-Viro model that accounts for the volume of 3-simplices was repro- 
duced correctly: as the geometrical considerations we used above tell us, we 
considered only the terms that are relevant for the volume of individual 3- 
simplices. However, there are other types of terms that we did not consider 
and that may be crucial to reproduce the corresponding state sum model 
correctly. Let us consider, for example, the terms of the type: 

Jtt(E Wi A E Wl A E W2 ), 

for some w\ 7^ wi- The result of such an integral is ill-defined. Indeed, there 
is an indeterminacy of the type • 00, where comes from du\ A du\ (see 
(f27|)), and 00 comes from the square of the 5- function 5 2 {u\). These terms 
are proportional to 

Tr ( X Wl X Wl X W2 ) . 

Other terms that may arise are 

jTr(E w AE w AE w ), 

for some w. They arise as the result of the indeterminacy 2 • oo 2 , where 
2 comes from du A du A du (see (p7|)), and oo 2 comes from the cube of the 
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5- function 5 3 (u). These terms are proportional to 

T^{X W X W X W ). 

There is no obvious reason to set these two types of terms to zero. In 
fact, as one can see, these are exactly the types of terms that are needed to 
account for the other terms appearing in the first A-order in the Turaev-Viro 
model. Although the structure of terms appearing this way is clear, at the 
present state of the development of our approach, no geometrical arguments 
is available to fix the ambiguity in the coefficients in front of such terms. 
However, the precise agreement of the terms for which such arguments do 
exist gives hope that, ones the above ambiguities are resolved, we will have 
an exact agreement between the models obtained using our procedure and 
the usual state sum models. 



5.4 4D BF theory with cosmological term 

In the case of 4D BF theory with cosmological term the interaction i<3? is 
given by: 

i$(E) = ~ f tr(E A E). (107) 
2 Jm 

Our general prescription is to evaluate this interaction term on the distri- 
bution (^) and find a polynomial function &(X). Thus, one has to evaluate 
integrals 

/ Tr{E w AE w >) (108) 
Jm 

with E w given by (^G|). The integral ( |108| ) is non-zero only if wedges w, w' 
intersect. However, similarly to the case of 3D, the result when the wedges 
intersect is ill-defined because of the distributional nature of E w . Thus, 
again some independent considerations have to be used to fix the ambiguity. 
We use a strategy similar to the one adopted in the case of 3D BF theory. 
We consider only the terms coming from wedges uu, w' intersecting at the 
center of a 4-simplex. As we shall see, these are the most important terms 
in the sense that they are responsible for the main terms in the first order 
in A of the Crane- Yetter model. The relevance of other types of terms 



will be emphasized below. To fix the ambiguity in (108) when w,w' are 
two wedges that intersect at the center of a 4-simplex h we use geometrical 
considerations. 

Let us consider the integral 

- - / Tr(E A E) (109) 
2 Jh 
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over the interior of a particular 4-simplex h. It is equal to the sum 

-\ £ / ME W AE'J. (110) 

Each of the integrals here is proportional to 

Tr(X w X w ,) = -2X l w X l w ,, (111) 
with the proportionality coefficient given by: 

- - / 8{u)8(y)8(u')5{v')du AdvA du A dv'. (112) 
2 Jh 

The later is equal to — (l/2)a • sign(u>, w'), where a is a numerical param- 
eter, whose value is not fixed due to the ambiguity referred to above, and 



sign(w,w') is the sign of the volume form in ( 112 ). Thus, ( |109| ) is equal to 

J2 a-sign(w,w')X % w X l w „ (H3) 

w,w'£h 

where the sum is taken over wedges w, w' inside h that span a 4-volume. 
There are exactly 30 terms summed over in (113). 

We will fix the parameter a relating the quantities X^ to the geometrical 
4-simplex in 1R 4 . First, let us note that when E in ( |109| ) is equal to the self- 
dual part of the wedge product of two copies of the frame field: 

E a b = + ^abi 

= tfcPi], (114) 

where 0^ is a frame (tetrad) field, then ( |109|) is equal to 

f l+ v /jl+ v ~a 

J h 2 ab 2 cdIJ 



1+yii +v* 2 a ^ c< ^ / X +V^^ X +V -abed 

h 2 2 



^ J h e U KL^J^ d L e abcd = ^V h , (115) 



where Vh is the volume of h with respect to the metric defined by 9^. Thus, 
we will fix a in such a way that ( |113j ) is equal to (3/2)V/j when X l w can be 
related to the quantities characterizing a geometrical 4-simplex. 

Recall that a geometrical 4-simplex in IR 4 is characterized (up to trans- 
lations) by four vectors: the four vectors pointing from one of the vertices 
to the other four ( for more details see the Appendix E). For each face of h 
one can also construct the so-called bivectors, which are given by the wedge 
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products of any two of the edges belonging to that face. The bivectors live 
in the second exterior power of IR 4 . One can take the self-dual part of a 
bivector to obtain an element of IR 3 . As is shown in the Appendix E (see 
( |166j )), the norm (obtained using the usual flat metric in K 3 ) of the self-dual 
part of each bivector is equal to the squared area of the corresponding face. 



There exists also a simple expression (167) for the volume of a 4-simplex in- 



volving only the self-dual parts of bivectors. Recall now that wedges w 6 h 
are in one-to-one correspondence with faces of h. Then, if has the in- 
terpretation of the self-dual part of the bivector corresponding to a face of 
h. 



A comparison of (167), (|115|) fixes the value of a to be a = 1/30-4. 



Thus, the first A-order term of the transition amplitude is given by: 

A S5 <»«> 

To compare this with (^) we need the relation between the grasping there 
and the operator in (|116|). The corresponding relation is given by 



l_6 5__ 

= 2SiJisiJi, {n } 

e' 

Using this correspondence we see that (|116| ) agree with the result (|SlJ). As 
in the 3D case, the matching includes the numerical coefficient and the sign 
of the expressions. 

Let us now discuss the role of the terms denoted by dots in (|8l|). Let us 
recall that those terms are determined by a framing of the (15j)-symbol used 
in the Crane- Yetter model. These terms are given by a sum of graspings of 
edges of the graph Th that share a vertex. Recall that vertices of Th are 
in one-to-one correspondence with the tetrahedra of h. Thus, the general 
structure of these terms is such that they can be grouped according to a 3- 
simplex (tetrahedron) to which they "belong" . Therefore, these terms can be 
thought of as the contribution to the 4- volume of the BF theory coming from 
individual tetrahedra of A. Thus, these terms are in certain sense analogous 
to the terms in 3D that carry 3- volume corresponding to the edges of A. 
As in 3D our approach did not reproduce correctly the terms corresponding 
to edges, in 4D we did not reproduce correctly the contribution to the 4- 
volume coming from tetrahedra. Similarly to the case of 3D, these terms may 
possibly be reproduced if one takes into account the types of intersections 
in (|108|) other than the ones considered here. 
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5.5 4D Euclidean gravity 



In order to apply our strategy to the self-dual Plebanski model let us consider 
the "interaction" term of the action (|43|): 

i${E) =~ij i>H A Ei - ^5 ij E k E^j . (118) 

Using the same procedure as in the case of the 4-dimensional BF theory with 
cosmological constant, this term, when evaluated on distributional E field, 
gives 

-*5>tf (119) 
h 

where denotes the center of the 4-simplex and 

^(X) = Y, sign(«;,«/)f^i/-^X^'fc) (120) 



w,w'(zh 



is a quadratic function in the wedge variables X w . In order to write (|119|) 
we have absorbed some unimportant numerical constants into the Lagrange 
multipliers tpij. The generating functional approach tells us that the spin 
foam transition amplitude of the self dual Plebanski model is given by Z(0) 
where 

Z(J)= ([[^(^e^^feWfil?).^/). (121) 
•* h 

The integration over ip 13 gives rise to delta functions. Thus, Z(J) can be 
characterized as the solution of 

^(4j)Z(J) = 0, (122) 
id J 

for all 4-simplices h, or, equivalently, as the solution of 

n ^4j ] ■ • • n ^4j )2{j)l j= ° = °- (123) 

Using the fact that Z(J) and, therefore, Z{J) are gauge invariant and the 
property that an su{2) symmetric traceless tensor fP J is totally characterized 
up to gauge transformation by its square we see that the preceding 

equation is equivalent to 

(^^) n (^7)-^)U=o = 0. (124) 
The solution of the later can be written as 

2(0) = hmn^e-^^^K^T) • Z(J)\j =0 (125) 
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In this form the result is very similar to the one given by the Reisenberger 
model. There are however several differences. The first one is the fact that 
the derivative operators appearing in our result are the commutative ones, 
while Reisenberger uses the non-commutative right and left invariant vector 
fields. The second difference is due to the presence of the factors of P(J) in 
the generating functional. In the next section we will emphasize significance 
of these discrepancies. 



6 Discussion 

In this section we summarize what has been achieved, discuss open questions 
of the approach and give some speculations as to possible directions of future 
research. 

As we have said in the introduction, our main aim was to understand 
whether the known spin foam (state sum) models, in particular the ones 
discussed in section |], could be obtained as a result of some systematic pro- 
cedure that starts from the corresponding classical action principle. We have 
proposed such a procedure, motivated by certain results of the loop approach 
to quantum gravity. Our main idea is two-fold. First, we proposed to include 
in the path integral certain distributional field configurations. In fact, in this 
paper we concentrated only on such distributional configurations. The dis- 
tributional field was taken to be the E field. As we discussed above, this is 
motivated by the distributional nature of the non-abelian "electric" field of 
the canonical approach. We put the distributional E field on a collection of 
two-dimensional polygons in spacetime, for which we used the collection of 
dual cells (faces) of a fixed triangulation of the spacetime manifold. Second, 
to calculate the path integral, we have employed the idea of the generating 
functional. This allows one to discuss different theories from the same point 
of view. As we have seen, the generating functional is the same for any the- 
ory in a given spacetime dimension, and transition amplitudes for different 
theories arise as different formal power series in the derivatives with respect 
to the current. Thus, this is quite reminiscent of the standard quantum field 
theory in Minkowski spacetime picture, where many different theories can 
arise from the same free theory as sums over different types of Feynman dia- 
grams. As we have seen, the role of the free theory is played in our approach 
by BF theory. Below we shall speculate more on the analogy between the 
diagrammatic expressions arising in our approach and Feynman diagrams. 

We have seen that our approach reproduces correctly the transition am- 
plitudes for BF theory in any dimension and the partition function for Yang- 
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Mills theory in two dimensions. We have also compared the transition am- 
plitudes for BF theories with cosmological term that arise in our approach 
with the ones given by the corresponding states sum models. We did this by 
comparing the first order terms in the decomposition of the corresponding 
amplitudes in the power series in A. In the case of 4D theory, we found 
agreement (apart from the terms that have to do with the framing), and 
in the 3D case only a partial agreement was reached. Note, however, that 
our approach does account for the most intricate part of the 3D and 4D 
transition amplitudes in the first A-order, reproducing correctly even the 
numerical coefficients. We have also seen, both in the cases of 3D and 4D, 
that the missing contributions can be interpreted as coming from singular 
contributions arising when one regularizes the volume operator. We expect 
that some new ideas will be necessary to resolve these regularization ambi- 
guities. Presumably, a study of the higher order terms in A may give some 
clues as to how these ambiguities must be resolved. 

Let us conclude with some speculations and directions for future re- 
search. An interesting open question that arises within our approach is as 
follows. We have seen that the state sum formulation of the topological 
BF theory can be obtained by using distributional fields associated with a 
2-dimensional cellular complex. Moreover, the residual regularized action 
on the 2-dimensional worldsheet is given by 2D BF- Yang-Mills theory. It is 
important to understand why such result is valid. One possibility would be 
to use the degenerate solutions found by Baez Q as some monopole contri- 
butions playing a major role in defining the measure of the field theory. It 
might be that the topological nature of the theory depends on whether these 
singular solutions dominate or not the partition function. 

In this paper we have discussed only correlation functions in the E fields. 
It is also possible to use our approach to compute the expectation value of 
the Wilson loop functionals. We hope to perform such calculations in the 
future. 

Let us now comment on the analogy between the diagrammatic expres- 
sions we have found and Feynman diagrams. In our approach the role of 
the "free" theory is played by the BF theory. Thus, the free transition am- 
plitudes in our case are just the BF transition amplitudes, given by certain 
sums of products of functions of spins labelling the triangulation. This is 
to be compared with the Feynman diagram technique free amplitudes which 
are just lines representing the free propagation of quanta. It is important 
that the BF amplitudes are independent of the triangulation used to com- 
pute them. They depend only on the properties of the initial and final states 
(and, of course, on the topology of the manifold that is contained between 
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the initial and final hypersurfaces). The more complicated, "interacting" 
amplitudes are given in our approach by inserting various chord diagrams 
into spin networks, with spin networks representing the free amplitudes. 
This is to be compared with how one takes into account interactions in the 
standard QFT's: graphically this is also represented by inserting certain 
types of vertices into the free diagrams. In our case, each grasping occurs 
within a particular simplex of the triangulation. One then has to take a sum 
over simplices. Similarly, in QFT's (in the coordinate representation) one 
puts an interaction vertex at a particular point and then integrates over the 
position of this point. Thus, it is tempting to draw an analogy between a 
simplex of a triangulation in our case and a point in QFT's. To summarize, 
spin networks play in our theory the role of the diagrams of free propagation, 
and the graspings by the chord diagrams are analogous to the insertion of 
interaction vertices into the diagrams of QFT's. 

Finally, let us discuss implications of the results we have obtained for 
the known spin foam models of quantum gravity in 4D. The most popular of 
these are the Reisenberger model ||, discussed in some details above, and the 
Barrett-Crane model H. The major problem with both of them is that they 
are rather invented then derived as a result of some systematic procedure 
starting from the Einstein-Hilbert action or any equivalent of it. However, as 
we have seen on the example of the Reisenberger model, and as one can check 
for the case of Barrett-Crane model, our approach to some extent reproduces 
the transition amplitudes given by these models. For example, our approach 
gives a clear interpretation of the so far mysterious evaluation of the spin 
networks that satisfy the constraint equations on the flat connection that is 
performed to obtain the amplitude. In our approach this evaluation has an 
obvious meaning of taking J = at the end of the calculation of a transition 
amplitude. We have also seen that the constraint equations of the known 
spin foam models have the meaning of differential equations that must be 
satisfied by the generating functional. However, although the main features 
of the both models were reproduced by our approach, some important details 
are different in the models we have obtained. Thus, our approach prescribes 
to use the commuting variational derivatives with respect to the current 
instead of non-commuting vector fields used by Reisenberger Q. Also, the 
presence of the function P(J), which plays a crucial role in, for example, 
2D Yang-Mills theory, is not accounted for by the Reisenberger and Barrett- 
Crane models. 

One more important lesson that one can learn from our results is as fol- 
lows. On examples of 3D and 4D BF theories with cosmological term we 
have seen that the simplest regularization of the "interaction" (volume) term 
of the action that we have discussed above does not account for all terms 
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one gets from the corresponding "exact" models. However, the Reisenberger 
model, for example, uses exactly the same way of regularizing the constraints 
of the theory, for, using the terminology of our approach, it takes into ac- 
count only the simplest types of intersections of wedges when regularizing 
the constraint part of the action. Our experience with the BF theories tells 
us that this cannot be correct. Thus, it could be that the problems with 
the Reisenberger model (for example, no interesting solutions of his con- 
straint equations is known) arise simply because this model does not take 
into a proper account other terms that arise in the regularization of the 
constraints. Similar is true for the Barrett-Crane model. 

The original motivation for our research program was to find a systematic 
procedure that gives one a way to derive spin foam models of various theories 
starting from the corresponding action functionals. However, the procedure 
proposed in this paper is still far from being able to claim to have achieved 
this goal. Yet, by studying the procedure proposed, we have learned a lot 
about existing models and obtained an important insight as to their deep 
internal structure. Our results also led us to the conclusion that the existing 
models of quantum gravity in 4D are too oversimplified. In spite of these 
rather negative conclusions, we believe that some of the insights we have 
obtained in this paper will become a part of the final, possibly much more 
elaborate picture. 

Acknowledgements: We are grateful to John Barrett, John Baez, 
Matthias Blau and Jose Zapata for discussions. K. K. was supported in 
part by the NSF grant PHY95-14240, by Braddock fellowship from Penn 
State and by the Eberly research funds of Penn State, L.F was supported 
by the CNRS and a NATO grant. 

A Conventions and notations 

We use the following conventions for differential forms and integrals: 



where e ai '" ad is the Levi-Civita density, taking the values plus minus one 
in any coordinate system. Our convention for the curvature form is: F = 
dA + A A A. 

Throughout the paper the following symbols will stand for the following 



A = -.A ai ... a dx ai A ... A dx ap 



(126) 



pi 




(127) 
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elements of the triangulation 

e — for an edge 
/ — for a face 
t — for a tetrahedron 
h — for a 4-simplex 

We will also use the symbol e to denote edges of the dual complex (dual 
1-cells), and a to denote dual faces (dual 2-cells). 

All traces that we use in this paper are in the fundamental representation. 

Our SO (3) index conventions: 

X = ir i X i , J=^ T i J i . (128) 



"- iJk ■ (129> 



The quantity dim g (j) is the so-called quantum dimension of j dim g (j) 
[2j + l] q , where [n] q is the quantum number 



qn/2 _ g-n/2 

gl/2 - q -l/2 ( 130 ) 



having the property that [k] q = 0. 



B Summary of facts on SU(2) 

Here we give a short summary of some standard facts about the group SU(2). 
One can parameterize an element g of SU(2) by vectors Z from the Lie 
algebra. The corresponding relation is given by the exponentiation map: 

g = e z = e ^" V! /2 ) ( 131 ) 

where n l is a unit vector n % rii = 1, tp is a real positive parameter, and a % are 
the usual Pauli matrices: 

(o 1 ■a j ) = ie ijk a k + 5 ij . (132) 



Thus, ipn 1 is an element of IR 3 ~ su(2), and (|131| ) gives the exponentiation 
map. As one can easily check, 

= coe(V»/2)+»nVsin(^/2). (133) 
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Thus, to cover the whole SU(2) just ones, the parameter ip should takes 
values in the range: ip £ [0,47r]. 

The Haar measure on the group can be related to the usual Lebesgue 
measure in IR 3 by introducing a function P(Z) on the Lie algebra: 



i'/2 

Then P 2 (Z)dZ/32ir 2 gives the normalized Haar measure on SU(2) in terms 
of the Lebesgue measure dZ. 

The characters of irreducible representations are given by: 

s *w+im) t (135) 

sm(^/2) 

where j are half- integers (spins). 

The Fourier transform on SU(2) maps any function on the group into a 
function on the space dual to the Lie algebra. Let / be a coordinate on the 
space su(2)*. Then the Fourier transform is given by: 

Hf) •■= J ^-P(Z)e-*^<P(expZ). (136) 
The inverse Fourier transform is given by: 

<KexpZ) = J2^^- ) Jd^ lf(Z) kf), (137) 
j K ) i 



where the integrals are taken over the co-adjoint orbits - spheres of radius 
j + 1/2, and the measure d£lf on each orbit is normalized so that 



(Kl f = dwij =2j + l. (138) 
A particular case of (137) is the following simple formula for characters ( |135| ): 

x ' {eZ) = p(z)l dnfelf{Z) - (139) 



This is the famous Kirillov formula [27]. 



C Haar measure, intertwiners and spin networks 



Let us denote by Vi the spin i representation of SU(2) and by V* the dual 
representation, which in the case of SU(2) is isomorphic to V{. Let us denote 



L. FREIDEL, K. KRASNOV 



1237 



by 6j the corresponding isomorphism, and by Ri{g) the representation of the 
group element g in Vi- 

Throughout the paper we use the normalized Haar measure 

J dg = 1. (140) 

Let K be an intertwiner between the trivial representation and a rep- 
resentation V, that is K G i/omG(m, V). We will denote by K its dual 
X G HomG(V,lR*). Intertwiners are the basic building blocks of the so- 
called spin networks. Another usage of the intertwiners is to express the re- 
sult of integration of a product of group elements. Let us denote by K t ^'-' tri 
an orthonormal basis of Home (1R, Vi x ® ■ ■ ■ <8> Vi n )- The intertwiners K£ ''" ,! " 
have the property that 

K a Kp = 6 a> p, 

where the product between K and K is defined by the duality bracket V ® 
V* — > 1R. Then the integral of n group elements is given by 

f dgR tl (g)®...®R tn (g)=J2K 1 - i "Kl^ (141) 

where R%{g) is considered as an element of Vi <8> V* . To integrate a product 
where both g and g~ l appear one has to use the duality relation Ri(g) = 
Ri{9~ 1 ) = e iRi(9) e i- F° r instance 

f — 1 - 

J dg (Ri)mn {Rj)m'n' = ~j . ~ &ij ^mm' ^nn' ■ (^^2) 

This equality can be conveniently expressed graphically if one represents 
a matrix element (Ri) mn by a line labelled by i with the two open ends 
corresponding each to one of the indices m, n. Let us symbolically denote 
the integration over the group by a circle going around the lines representing 
the group elements. Then (|142j) takes the following form: 



<&> = ^ s -> ^ (143> 

To find a graphical representation of the result of the integrals involving 
more than two matrix elements we introduce a trivalent vertex - analog of 
the Clebsch-Gordan symbol - normalized in a special way. But first, let us 
define the so-called spin networks. 
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An SU(2) spin network functional is denned by the following data: (i) 
an oriented graph T; (ii) a map j from the set of edges of T to the set of 
irreducible representations of SU(2); (iii) a map i from the set of vertices 
V to the set of intertwiners, which assigns to each vertex v an intertwiner 
from the tensor product of representations labelling the incoming edges to 
the tensor product of representations labelling the outgoing edges. A spin 
network is a function on a number of copies of SU(2), more precisely, on 
(SU(2)) E , where E is the number of edges in T. To find the value of this 
function one takes for each edge the matrix element of the group element 
on that edge in the representation that labels the edge, and contracts the 
matrix elements corresponding to all edges using intertwiners at vertices. 
The function obtained this way is gauge invariant. 

In this paper we will use the normalized trivalent vertex, defined in such 
a way that the ^-symbol constructed from 3 spins is always equal to unity: 



where the operation (-)o denotes the evaluation of the corresponding spin 
network on all group elements equal to the unity in the group. 

With this normalization of the trivalent vertex intertwiner, one can show 
that the following relation holds: 



Using this relation one can find the result of the integral of a product of 3 
and 4 group elements. The corresponding formulas are: 




(144) 




(145) 



k 




(146) 



c: 




i 



k 



1 




V ■ V 



(147) 



A 




m 



m 



Also, we use the same trivalent vertex to define the normalized (6j),(15j)- 
symbols used in the body of the paper. Both symbols are given by evalu- 
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ations of the corresponding spin networks, where the intertwiner used in a 
trivalent vertex is always the one normalized as in ( |144|) . Thus, we get: 



(6j) 




(148) 



(15j) 



where the resolution of the 4-valent vertex in ( |149| ) is given by: 

5 j 



(149) 




(150) 



D Quantum 6j symbol 



Let X be a one dimensional oriented compact manifold, or more generally 
an oriented graph. A chord diagram (usually referred as Chinese character 
chord diagram) with support X is the union D = DU X where D is a graph 
with univalent and trivalent vertices, together with a cyclic orientation of 
trivalent vertices and such that univalent vertices lie in X. Trivalent vertices 
are referred to as internal vertices, and the degree of D denoted by d°(D) 
is half the number of vertices of the graph D. Let A n be the Z module 
freely generated by chord diagrams of degree n. We define the Z module 
of Vassiliev diagrams of degree n denoted by A n as being the quotient of 
A n by the relations (STU, IHX, AS) shown in Figure ||. In the figures we 
always represent the support X with bold lines and the graph D with dashed 
chords. 

Given a Lie group G, and a coloring C of the graph it is possible to 
define an evaluation of chord diagrams that is usually called weight system 
and denoted by lug,c- Here C, the coloring of the graph X, is a map which 
assigns a representation of G to any edge of the graph and an intertwiner 
to any vertex of the graph X. Given a chord diagram D with support X 
and a coloring C of A, we define u>g,c{D) as shown in Figure ^. Here X a 
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± \ w — I 



) ^ - /- - - N = 1 

/ \ / \ A 

/ \ / \ 




= 



Figure 8: 



denotes a basis of the Lie algebra of G, f denotes the structure constants 
associated with each internal vertices and t = t ab X a (g) Xb is the quadratic 
casimir t = t ab X a (g) Xf, associated to each chord. 




b c 



Figure 9: 

The space of chord diagrams has been used to define the now famous 
universal invariant of framed link Z p8| , which assigns to any link L in i? 3 
a formal power series in h Z(L) = ^2 n h n Z n (L), with Z n (L) G A n . We are 
not going into the details of this construction. However we must say that 
this construction uses three main building blocks: The braiding R G A(I 2 ), 
the associator (p 6 ^I 3 ), which should satisfy the pentagon and hexagon 
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identity, and the value of the un-knot v G A(I), where / denotes the unit 



interval (|15T|) J29 



12 



R 



12 



+ 



+ ••• 



(151) 



/ / 



+ • 



(152) 



(153) 



One important property of the universal invariant is that log,c{Z{L)) is 
equal to the Reshetikhin-Turaev evaluation of the link L associated with 
the quantum group U g (G) 3C, 32], where q = expift and the casimir is 



normalized such that the norm of long roots is 2. 

It has been shown recently |34| ], that the universal invariant Z can be 
extended to an invariant of oriented, framed trivalent graphs embedded into 
i? 3 . The evaluation of the tetrahedral graph with the universal invariant is 



given by equation ( 154 ). 





(154) 



If we evaluate the result of Fig. |154| using the weight system uj su (2),c we 



get the normalized Gj symbol expanded in terms of fi (Reshetikhin-Turaev 
evaluation of the tetrahedral graph for \J e m (su(2))). The term proportional 
to tt° is just the classical 6j symbol. The next term, proportional to h 2 , is 
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Figure 10: A 4-simplex is characterized up to translations by four vectors. 



given by : 






(155) 



Here we have used the fact that, in the case of su{2) and for the normalization 
of the casimir given by the trace in the fundamental representation, the 
identity = 4 is satisfied. This term can be written in a symmetric form 
as: 



tet 



(i) 



-<- £< 

16 v 24 V„ 




)), (156) 



where the first sum is over all triples of distinct edges of the tetrahedra and 
the second sum is over all edges of the tetrahedra. 



E Geometric 4-simplex 



In this Appendix we list some facts about the geometry of a 4-simplex in 
IR 4 . The geometrical considerations used in ( |5.4| ) are based on some of these 
facts. 

A 4-simplex in IR 4 is characterized (up to translations) by four vectors. 
These, for instance, can be vectors pointing from one of the vertices, which 
we will denote by (0), to the other four vertices (l)-(4). See Fig. 10. Let us 
denote these vectors by e„. Here I = 1 ... 4 is an index for a vector in IR 4 , 
and a = 1 ... 4 indicates a vertex at which the vector is directed. Thus, e[ 
is a vector pointing from vertex (0) to the vertex (1). 
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Instead of vectors it is sometimes more convenient to use the so-called 
bivectors. A bivector E IJ is an element of the second exterior power of IR 4 . 
In particular, it can be obtained by wedging two vectors. Thus, bivectors 
that characterize a 4-simplex are obtained as E^ = e^e^ . Here the brack- 
ets denote the operation of antisymmetrization. It is not hard to see that 
bivectors E^ are in one-to-one correspondence with faces of the 4-simplex 
h. For example, the bivector E\^ corresponds to the face whose vertices are 
(0),(1),(2). The norm of each bivector is proportional to the squared area of 
the corresponding face 

E'jEabjj = 2A 2 ab , (157) 

where A a b is the area of the face (0),(a),(b) and no summation over a, b is 
assumed. The volume of h can be obtained by wedging two bivectors that 
correspond to faces that do not share an edge. For example, with our choice 
of orientation of M, the volume is given by 

V h = ^e IJKL E{iEg L . (158) 



It is sometimes convenient to introduce bivectors corresponding to all 10 
faces of h. So far we have introduced 6 bivectors #£/ corresponding to 6 
faces of h. Bivectors that correspond to other 4 faces can be obtained as 
linear combinations of these 6 bivectors. When working with all 10 bivectors, 
it is convenient to label bivectors by 3 different indices instead of just two. 
We will employ for this purpose small Greek letters. Thus, bivectors are 
denoted by E 1 ^ , a, (3, 7 = ... 4. These bivectors are defined by 

E*&y = i>2y> (159) 

where e£g is a vector that points from the vertex a to vertex j3. The norm 
of all bivectors ( |159|) is equal to the twice of the squared area of the corre- 
sponding face, as in Q157 ). As we have said above, only 6 of 10 bivectors 



(|15S| ) are independent. Thus, there are certain relations between them. One 
can write one such relation for each tetrahedron of h. One gets 5 relations 
only 4 of which are independent. With our definition (|159| ) these relations 



arc 



#012 + #023 — #013 — #123 = 0, (160) 
#013 + #034 — #014 — #134 = 0, 
#024 + #234 — #023 _ #034 = 0, 
#014 + #124 — #012 — #024 = 0, 
#123 + #134 — #124 — #234 = 0, 
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where we have suppressed the indices /, J for brevity. The volume of h can be 
expressed as a wedge product of any two of the bivectors (|159|) corresponding 
to faces that do not share an edge. This can be written as 

sign(/, f')V h = ±e IJKL E{f) IJ E{fi) K \ (161) 

where we have introduced a notation E(f ) IJ for a bivector that corresponds 
to face /, and sign(/, /') is the sign of the right-hand-side in (|161| ). The 
expression ( |161| ) gives the volume of h for any two pairs of faces /, /' that 
do not share an edge. One can use the expression (|158| ) and the "closure" 
relations ( |160| ) to work out the correct sign of any of such formula for the 
volume. 

Any bivector can naturally be split in its self-dual and anti-self-dual 
parts: 

E a f3y = + E a f3 y + E a p T (162) 
The self-dual and anti-self-dual parts are given correspondingly by 

+ E a p 1 = -(Ea/3^ + *E a py), (163) 

EoiP-f = -^{Eap-i — *E a p 1 ), 

where the Hodge star duality operation is defined as 

* E Li-r = 2 e KL E a^c (I 64 ) 

Since the space of self-dual (and anti-self-dual) bivectors in 1R 4 is three- 
dimensional, we can introduce a new set of indices to label them. Thus, as 
the index for self-dual (anti-self-dual) bivector we will use lower case Latin 
letters from the middle of the alphabet: k, . . . = 1, 2, 3. The norm of any 
self-dual (anti-self-dual) bivector calculated by contracting indices /, J will 
be the same as the norm calculated by contracting the single index v. 

+ E i+ E, t = + E ij+ Ejj, (165) 

where we have suppressed the indices a, /3, — 

Not any bivector in IR 4 is simple, that is, not any bivector is a wedge 
product of two vectors. The necessary and sufficient requirement of sim- 
plicity is that the norm of the self-dual part is equal to the norm of the 
anti-self-dual part of the bivector. Thus, using ( |157| ), we can conclude that 
when a bivector is simple, its self-dual part norm is equal to the squared 
area of the corresponding face. 

+ Ea / 3 J + E a p li = A 2 af}l . (166) 
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There exists an expression for the volume of h that involves only the self-dual 
parts of bivectors ( |159|) : 



Vfc = ~£sign(/,/0 + £W^(A (167) 



where + E l (f) is the self-dual part of the bivector ( |159|) corresponding to 
the face /, the sum is taken over all pairs /, /' of faces that do not share 



an edge, and sign(/, /') is the function introduced above by equation ( 161 ). 



The factor of 1/30 in (167) appears to cancel the factor that comes from the 
sum over 30 terms /, /'. 
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